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We treat a relativistically moving particle interacting with 
a quantum field Irom an open system viewpoint of quantum 
field theory by the method of influence functionals or closed- 
time-path coarse-grained effective actions. The particle trajec- 
tory is not prescribed but is determined by the backreaction of 
the quantum field in a self-consistent way. Coarse-graining the 
quantum field imparts stochastic behavior in the particle trajec- 
tory. The formalism is set up here as a precursor to a first princi- 
ples derivation of the Abraham-Lorentz-Dirac (ALD) equation 
from quantum field theory as the correct equation of motion 
valid in the semiclassical limit. This approach also discerns clas- 
sical radiation reaction from quantum dissipation in the motion 
of a charged particle; only the latter is related to vacuum fluctu- 
ations in the quantum field by a fluctuation-dissipation relation, 
which we show to exist for nonequilibrim processes under this 
type of nonlinear coupling. This formalism leads naturally to 
a set of Langevin equations associated with a generalized ALD 
equation. These multiparticle stochastic differential equations 
feature local dissipation (for massless quantum fields), multi- 
plicative noise, and nonlocal particle-particle correlations, in- 
terrelated in ways characteristic of nonlinear theories, through 
generalized fluctuation-dissipation relations. 



I. INTRODUCTION 



X 



A. Particles and fields 

Charged particles moving in a quantum field is an old 
topic in electromagnetic radiation theory, plasma physics 
and quantum/atom optics. Interesting features include the 
relation between quantum fluctuations and radiation re- 
action 0, and, for the case of uniformly accelerated de- 
tectors, thermal radiance in the detector (known as the 
Unruh effect @]). For relativistic particles there are also 
bremsstrahlung, synchrotron radiation, and pair creation. 
The motion of congruences of charged particles shows up 
in particle beam, plasma and nuclear (e.g. heavy-ion colli- 
sion) physics. Their collective or transport properties must 
be treated by additional statistical mechanical considera- 
tions beyond those applied to single particles. 

Theoretically, since particle-field interaction is in princi- 
ple described by quantum field theory, one may get the im- 
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pression that ordinary field theoretic methods are both nec- 
essary and sufficient for the study of this problem. Quan- 
tum field theory, in the way we usually learn it, is customar- 
ily formulated in the context of, and with special emphasis 
on, how to answer questions posed in particle physics. For 
example, the S matrix is constructed for calculating scat- 
tering amplitudes, a perturbation expansion (e.g. Feynman 
diagram techniques) is used for theories with a small cou- 
pling constant (e.g., QED), and even the effective action 
is usually based on in-out (Schwinger-DeWitt) boundary 
conditions (the vacuum persistence amplitude). But, when 
one wants to find the evolution equation of a particle one 
needs an in-in boundary condition and an initial value for- 
mulation of quantum field theory ||]4) . This is the starting 
point for constructing theories describing the nonequilib- 
rium dynamics of many-particle systems |g| where both the 
statistical mechanics depicting the collective behavior of a 
congruence of particles as well as the statistical mechanical 
properties of interacting quantum field theory need to be 
included in our consideration 1^,0. 

The relationship between the contrasting paradigms of 
particles and fields is at the heart of our inquiry. The con- 
cept of a particle is very different from that of a quan- 
tum field: a particle moving in real space cannot be com- 
pletely described by, say, a number representation in Fock 
space, which underlies the second quantization formula- 
tion of canonical field theory. Background field decompo- 
sition into a classical field and its quantum fluctuations 
is a useful scheme, where one can follow the development 
of the classical background field with influences from the 
quantum fluctuations. While these methods have been ap- 
plied extensively to the semiclassical evolution of quantum 
fields, the classical picture of particle motion is still quite 
remote. The Feynman path integral formulation (for par- 
ticles worldlines as opposed to fields) makes it easier to 
introduce classical particle trajectories as the path giving 
the extremal contribution to the quantum action, the sta- 
tionary phase defining one condition of classicality. One 
can then add on radiative corrections by carrying out the 
loop expansion. This approach puts a natural emphasis 
on particle trajectories rather than scattering amplitudes 
between momentum eigenstates. But how does one go fur- 
ther to incorporate or explain stochastic behavior of the 
particle? We shall see that when the coupling between 
the particle and the field is non-negligible, backreaction of 
the field (with all its activities such as vacuum polariza- 
tion, pair creation, etc.) on the particle (beyond simple 
radiative corrections) imparts stochastic components to its 
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classical trajectory. 

At a deeper level, which we will tend to in the next series 
of papers, even a simple quest to understand the detailed 
behavior of quantum relativistic trajectories (i.e. world- 
lines) is a complex challenge, with problems connected with 
many central issues of quantum physics. For example, un- 
derstanding how a particle moves through time (or more 
accurately, how a property that can be recognized as time 
emerges for a single particle) provides technical and concep- 
tual insights into the larger questions concerning the nature 
of time in general. Here we focus on the regime where the 
notion of trajectory and particle identity is well-defined- 
this implies that both particle creation/annihilation and 
quantum exchange statistics (e.g. Fermi exclusion or Bose 
condensation effects) are of negligible importance. 

Is there a unified framework to account for all these as- 
pects of the problem? This is what we strive to develop 
in this series of papers. We see from the above cursory 
queries that to address this deceptively simple problem in 
its full complexity one cannot simply apply the standard 
textbook recipes of quantum field theory. Even a "sim- 
ple" problem such as particles moving in a quantum field, 
when full backreaction is mandated, involves new ways of 
conceptualization and formulation. In terms of new con- 
ceptualization, it requires an understanding of the relation 
of quantum, stochastic, classical (as ingrained in the pro- 
cess of decoherence) behaviors, and of the manifestation 
of statistical mechanical properties of particles and fields, 
in addressing where dissipation and fluctuations arise and 
how the correlations of the quantum field enter. In terms 
of new formulation, it involves the adaptation of the open- 
systems concepts and techniques in nonequilibrium statis- 
tical mechanics, and an initial value (in-in) path integral 
formulation of quantum field (and quantized worldline) the- 
ory for deriving the evolution equations. We will discuss 
these two aspects to highlight the basic issues involved in 
this problem. 

In this investigation we take a microscopic view, using 
quantum field theory as the starting point. This is in vari- 
ance to, say, starting at the stochastic level with a phe- 
nomenological noise term often-times put in by hand in 
the Langevin equation. We want to give a first-principles 
derivation of moving particles interacting with a quantum 
field from an open-systems perspective. A consequence of 
coarse-graining the environment (quantum field) is the ap- 
pearance of noise which is instrumental to the decoherence 
of the system and the emergence of a classical particle pic- 
ture. At the semiclassical level, where a classical particle is 
treated self-consistently with backreaction from the quan- 
tum field, an equation of motion for the mean coordinates 
of the particle trajectory is obtained. This is identical in 
form to the classical equations of motion for the particle 
since higher-order quantum effects, which arise when there 
are nonlinear interactions, are suppressed by decoherence. 

Backreaction of radiation emitted by the particle on the 
particle itself is called radiation reaction. For the spe- 
cial case of uniform acceleration it is equal to zero. This 



well-known, but at first-sight surprising, result is consis- 
tent owing to the interplay of the so-called acceleration 
field and radiation field ||. Radiation reaction (RR) is 
often regarded as balanced by vacuum fluctuations (VF) 
via a fluctuation dissipation relation (FDR). This leads to 
a common misconception: RR exists already at the clas- 
sical level, whereas VF is of quantum nature. There is, 
as we shall see in this paper, nonetheless a FDR at work 
balancing quantum dissipation (the part which is over and 
above the classical radiation reaction) and vacuum fluctu- 
ations. But it first appears only at the stochastic-level, 
when self consistent backreaction of the fluctuations in the 
quantum field is included in our consideration. In addition 
to providing a source for decoherence in the quantum sys- 
tem making it possible to give a classical description such 
as particle trajectories, fluctuations in the quantum field 
are also responsible for the added dissipation (beyond the 
classical RR), and a stochastic component in the particle 
trajectory (beyond the mean). Their balance is embodied 
in a set of generalized fluctuation-dissipation relations, the 
precise conditions for their existence we will demonstrate 
in this paper. 

Ultimately a satisfactory description of the particle-field 
system would have to come from a full quantum theory 
treatment. One needs to stipulate and demonstrate clearly 
what successive approximations one introduces to some 
larger interacting quantum system will enable us to begin 
to see the precursor or progenitor of the particle, and the 
residual or background quantum field, depicted so simply 
in the end as the particle-field system in the conventional 
field theory idiom. The imprint of the particle's motion is 
left in the altered correlations of the quantum field || . 

Later we shall see that while the problem of the nonequi- 
librium quantum dynamics of relativistic particles (or 
fields) is important in its own right; the rich physics con- 
tained is also closely connected with other problems of fun- 
damental significance such as Quantum Gravity and String 
theory, including semiclassical gravity and quantum fields 
in black hole and early universe spacetimes. The remain- 
der of this introduction summarizes the main ideas in this 
work and places it in the diverse range of prior research 
on this subject matter. We try to present a viewpoint and 
method which may provide a comprehensive and unified 
account while pointing out specific areas requiring further 
attention. 



B. Quantum, Stochastic and Semiclassical Regimes 

Our treatment in this first series emphasizes the semiclas- 
sical level and stochastic regimes. Our approach is sharply 
distinguished from the case where the motion of the particle 
is prescribed, i.e., a stipulated trajectory, that is character- 
istic of many treatments of 'particle-detectors' in quantum 
fields. Here, the trajectory is determined by the quantum 
field in a self-consistent manner. In the former case, there 
is a tacit presence of an agent which supplies the energy to 
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keep the trajectory on a prescribed course, its effect show- 
ing up in the radiation given off by a detector/particle un- 
dergoing, say, acceleration. The field configuration would 
have to adjust to this prescribed particle motion accord- 
ingly, without affecting the particle motion. In the latter 
case, the only source of energy sustaining the particle's mo- 
tion arises from the quantum field, and both the field and 
particle adjust to each other in a self-consistent manner. 
The former case has been treated in detail by Raval, Hu, 
and Anglin (RHA) |9| (see references therein for prior work) 
using the concept of quantum open-systems and the influ- 
ence functional technique. We investigate the second class 
of problems now with the same open-system methodology. 

A closed quantum system can be partitioned into several 
subsystems according to the relevant physical scales. If one 
is interested in the details of one such subsystem, call it the 
distinguished (relevant) system, and decides to ignore cer- 
tain details of the other subsystems, comprising the envi- 
ronment, the distin guis hed subsystem is thereby rendered 
as an open-system ]Kj. The overall effect of the coarse- 
grained environment on the open-system can be captured 
by the influence functional technique of Feynman and Ver- 
non 1 11 1, or the closely related closed-time-path effective 
action method of Schwinger and Keldysh gj. These are 
initial value formulations. For the model of particle- field 
interactions we study, this approach yields an exact, nonlo- 
cal, coarse-grained effective action (CGEA) for the particle 
motion p2| . The CGEA may be used to treat the com- 
plete quantum dynamics of interacting particles. However, 
only when the particle trajectory becomes largely well- 
defined (with some degree of stochasticity caused by noise) 
as a result of effective dccoherence due to interactions with 
the field can the CGEA be meaningfully transcribed into 
a stochastic effective action, describing stochastic particle 
motion §JlJ]JlJ]. 

The effect of the environment (the coarse-grained sub- 
systems) on the system (the distinguished subsystem) is 
known as backreaction. One form of classical backreaction 
in the context of a moving charged particle in a quantum 
field is radiation reaction which exerts a damping effect 
on the particle motion. As remarked before, it should not 
be mistaken to be balanced directly by an FDR with vac- 
uum fluctuations of the quantum field. The latter induce 
quantum dissipation and are instrumental to decohering 
the classical particle. Let us now further examine the role 
of dissipation, fluctuations, noise, decoherence |||7]] and the 
relation between quantum, stochastic and classical behav- 
ior @0. 

Decoherence or dephasing refers to the loss of phase co- 
herence in the quantum open system arising from the inter- 
action of the systems with the environment. Effective deco- 
herence brings about the emergence of classical behavior in 
the system which generally carries also stochastic features. 
Under certain conditions quantum fluctuations in the envi- 
ronment act effectively as a classical stochastic source, or 
noise. While noise in the environment is instrumental in de- 
coherence, decoherence is a necessary condition for the ap- 



pearance of a classical trajectory. In this emergent picture 
of the quantum to classical transition, there is always some 
degree of resultant stochasticity in the system dynamics 
J13[ . When sufficiently coarse-grained descriptions of the 
microscopic degrees of freedom are considered, nearly com- 
plete decoherence leads to negligible noise, and the classi- 
cal description of the world is complete. Moving back from 
classicality towards a description of more finely-grained his- 
tories^], the quantum effects suppressed by decoherence on 
macroscopic scales reveal themselves in the emergence of 
stochasticity. In this realm, decoherence, noise, dissipation, 
particle creation, and backreaction are seen as aspects of 
the same basic quantum-open-system processes [|l3 14 . If 
'minimal' additional smearing is not adequate to decohere 
the particle trajectories, the stochastic limit is not physi- 
cally meaningful because the quantum interference between 
particle histories continues to play a significant role. Even 
without a field's presence as an environment, the particle's 
own quantum fluctuations (represented by its higher order 
correlation functions) may be treated as an effective envi- 
ronment for the lower-order particle correlation functions 
(particularly the mean trajectories), so that a stochastic 
regime may still be realized when appropriately coarse- 
grained histories are considered |gJ^,Q. When averaged 
descriptions are considered such that the coarse-grained set 
of histories has substantial inertia, the quantum-fluctuation 
induced noise is negligible, and one moves from the stochas- 
tic to semiclassical domains. Note that there is now grow- 
ing recognition that the transition from quantum to clas- 
sical behavior |L|[l5| is characterized by a rich diversity 
of scales and phenomena, and it should therefore not be 
thought of as a single transition, but rather a succession of 
regimes that may, or may not, be well-separated in practice 

The view of the emergence of semiclassical solutions as 
decoherent histories also suggests a new way to look at 
the radiation-reaction problem for charged particles. The 
classical equations of motion with backreaction are known 
as the Abraham-Lorentz-Dirac (ALD) equations [Q. The 
solutions to the ALD equations have prompted a long his- 
tory of controversy due to such puzzling features as pre- 
accelerations, runaways, non-uniqueness of solutions, and 
the need for non-Newtonian initial data ||. It has long 
been felt that the resolution of these problems must lie in 
quantum theory But, this still leaves open the questions of 
when, if ever, the Abraham-Lorentz-Dirac equation appro- 
priately characterizes the classical limit of particle backre- 
action; how the classical limit emerges; and what imprints 
the correlations of the quantum field environment leave. 
We show that these questions, and the traditional para- 
doxes, both technical and conceptual, can be resolved in 



1 The finest-grained histories are just the skeletonized paths in 
the path integral formulation. Coarse-graining is achieved by 
integrating over subsets of these paths. 
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the context of the initial value quantum open system ap- 
proach. 

Indeed, since every possible fine-grained history is in- 
cluded in the path integral for the quantum evolution, 
there is no a priori reason to reject fine-grained runaway 
solutions as unphysical, nor is there any sense in which 
a particular fine-grained history pre-accelerates, since the 
fine-grained paths that appear in the path integral aren't 
causally determined by earlier events in any case. From 
this perspective, the appropriate questions to address are, 
which 'quasi-classical' coarse-grained solutions decohere 
more readily. Certainly it would be strange if runaway (or 
pre-accelerating) decoherent coarse-grained histories oc- 
curred with any appreciable associated probability. So the 
question that should be asked, in the context of how classi- 
cal solutions arise from the quantum realm, is whether de- 
cohered particle- histories are 1) solutions to the ALD equa- 
tion, 2) unique and runaway free, and 3) without preac- 
celeration on scales larger than the coarse-graining scale 
at which the trajectories decohere. We shall see that, in 
fact, the semiclassical solutions do satisfy these criteria, 
and therefore, one is entirely justified in using the ALD 
equation in the classical limit. One also sees that because 
this approach describes coarse-grained histories, the ALD 
equation as an approximation fails in the finest-grained 
quantum limit. 

Therefore, a fundamental understanding of the stochas- 
tic and semiclassical limits must be set in the context of 
the full quantum theory of particles and fields. In the sec- 
ond scries of papers, we further develop the worldline plus 
field framework to describe fully quantized relativistic par- 
ticles in motion through spacetime interacting with quan- 
tum fields and highlight its semiclassical, stochastic and 
quantum features. This framework for understanding rela- 
tivistic systems is important because, in the noncquilibrium 
dynamics of real particles, the localized nature of the par- 
ticle state is a prominent characteristic of the semiclassical 
limit, and this fact is not most naturally described by the 
usual perturbative field theory in momentum space. Fur- 
thermore, interaction, correlation, measurement and deco- 
herence invariably take place in both space and time. The 
incorrect treatment of "measurement" as an instantaneous 
occurrence leads to fundamental inconsistencies in the de- 
scription of physical initial states which continues to be a 
substantial obstacle to a deeper understanding of decoher- 
ence and nonlocal correlation. For this reason, addressing 
the stochastic (and quantum) limits in a relativistic frame- 
work treating space and time covariantly is a vital element 
in fully understanding the stochastic-semiclassical limit. 
By developing a new framework for relativistic particle- 
field quantum dynamics we will show how this alternative 
approach provides powerful tools for addressing both con- 
ventional issues, and for exploring new questions. 



C. Coarse-graining the particle 

Our problem, as well as many others from quantum and 
atom optics, provides a good example of where a quantum 
field (e.g. photons) acts as an environment in its influ- 
ence on an atom or electron system. There are, of course, 
physical contexts for which it is more appropriate to coarse- 
grain the particle degrees of freedom, whereby one obtains 
a coarse-grained effective action for the field. In this com- 
plementary view, matter plays the role of an environment 
for the field as a system. When both particle and field 
coarse-grainings converge to mutually decoherent sets of 
particle and field histories, one recovers the classical limit 
119 . These regimes are illustrated schematically in Figure 
1, where the field degrees of freedom are denoted by tp, and 
the particle degrees of freedom are denoted by z. 

D. Nonlinear Coupling 

Comparing with prior work on this subject matter using 
the same approach, the most closely related being that of 
Raval, Hu and Anglin ||], who derived the influence func- 
tional for n-detectors moving in a quantum field, the main 
distinct feature here is that the particle trajectory is not 
prescribed, but is a dynamical degree of freedom deter- 
mined self-consistently by the field. As such the coupling 
between the particle(s) and the field is of a more funda- 
mental nonlinear nature. One type of nonlinear coupling 
between the system and the environment in the quantum 
Brownian motion models pQ-p^] has been considered by 
Hu, Paz, and Zhang Their model contains nonlin- 

ear couplings that are polynomial in the field variables, 
whereas here we consider the opposite example, where the 
field variables are linear, but the system variables are not. 
This is the case for QED, which is our ultimate goal (to be 
discussed in Paper III Q of this series). 

Our treatment of nonlinear particle-field interactions is 
in sharp contrast with the existing body of work on the 
semiclassical and stochastic limit for linear QBM. We find 
that the requirement of self-consistency, together with the 
nonlinearity of the fundamental interactions, leads to a de- 
scription of the stochastic limit that is more unified and 
tightly constrained than that which has been described up 
to now. In fact, confronting nonlinear systems is a crucial 
next step in exploring the stochastic limit, and the deeper 
relationship between noise, decoherence, and the intricate 
evolution of correlations for systems with many degrees of 
freedom. Most studies of decoherence have invoked linear 
systems with an (externally) predetermined basis of states 
which are then shown to decohere. But a complete de- 
scription of decoherence must include the mechanism by 
which a decoherent basis of states is self-consistently (e.g. 
not externally) selected. Understanding this mechanism as 
a dynamical feature of nonlinear theories is an important 
problem that cannot be fully addressed within linear mod- 
els. 
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Another important difference between linear and nonlin- 
ear theories is in their respective equations of motion for 
correlation functions. For linear theories, these are just 
given by the classical equations of motion, and the quan- 
tum evolution may be reproduced by a probabilistic de- 
scription of initial conditions. An example of this is the 
Wigner function evolution for linear systems with initially 
positive definite distributions in phase spaces. The equa- 
tions of motion are then the classical ones, and negative 
values of the Wigner function never evolve in the future. 
One may then view the theory as dynamically equivalent to 
classical statistical dynamics (though of course the mean- 
ing of the Wigner functions is still quantum p^|). This 
fact is what makes the construction of a stochastic limit 
that bridges linear classical and quantum theories fairly 
straightforward. The same is not true of nonlinear the- 
ories. Their quantum equations of motion (even for the 
mean-field) are not equivalent to the classical ones. This 
implies that there are important qualitative differences in 
the quantum dynamics of linear versus nonlinear theories 
that are significant even in the semiclassical regime, and 
that will be missed by the study of linear systems alone. 

The important result in this first paper is the derivation 
of self-consistent Langevin equations for relativistic par- 
ticle motion in a quantum field starting from relativistic 
quantum mechanics. A stochastic description of the limit 
of nonlinear theories requires great care, and, in its most 
general form, will involve nonlinear stochastic differential 
equations where the statistics of the induced particle fluctu- 
ations are externally determined by the quantum statistics 
of the fields. 

Because there are inherent dangers in a phenomenologi- 
cal approach to nonlinear stochastic equations, it is crucial 
to work from first principles stressing the microscopic ori- 
gin of fluctuations. In this work, we treat an example of a 
nonlinear stochastic system in the regime where decoher- 
ence allows the expansion around the semiclassical solu- 
tion. The decoherent semiclassical solutions need not, and 
usually are not, equilibrium solutions. Hence, compared 
with the traditional derivation of the Langevin equations 
in the context of near-equilibrium linear response theory, 
this method is significantly more general. 

For Langevin equations the effect of the field is registered 
in the stochastic properties of the particles. The stochastic 
mean of the equations of motion corresponds to the mean- 
field semi-classical limit; but in addition, the symmetrized 
n-point correlation functions for the particles are given 
by the higher order stochastic moments. Therefore, the 
stochastic particles may be thought of as being "dressed" 
by the non-local statistics of the field. These equations 
impart stochastic features to the particle and field prop- 
erties beyond the semiclassical limit, leading towards the 
quantum domain. The stochastic equations of motion fea- 
turing nonlocal noise, causal particle-particle interactions, 
and nonlocal particle-particle correlations make the evolu- 
tion highly non-Markovian (i.e. history dependent). Only 
in the single particle semiclassical limit (with local dissipa- 



tion and white noise) is the evolution strictly Markovian. 
For well separated particles, in the high temperature limit, 
the Langevin equations may be approximated by Marko- 
vian dynamics; this is the regime in which the nonlocal role 
of the quantum vacuum has been essentially washed out. 
But as one moves more deeply into the stochastic regime 
(especially at low temperatures), the Markovian approxi- 
mation is no longer adequate. It is in this realm that our 
methods are markedly distinct from the Hamiltonian meth- 
ods more commonly used in the derivation of Markovian 
master equations. 

E. Prior work in relation to ours 

The treatment of a quantum field as a bath of har- 
monic oscillators has a long history. Many authors take 
a semi-phenomenological approach in adding a noise to the 
quantum equation of motion by hand to get the Quantum 
Langevin equations (QLE). This practice works reasonably 
well in the linear response regime for a system in equilib- 
rium, but can otherwise violate the fluctuation-dissipation 
relation and bring about pathologies (such as a-causal evo- 
lution). There are many ways to derive the QLE describing 
the dissipative dynamics of a quantum system in contact 
with a quantum environment, such as modified (Heisenberg 
picture) Hamilton equations of motion or path integrals. 
Caldeira and Leggett's revival of the Feynman- Vernon in- 
fluence functional in the study of quantum Brownian mo- 
tion (QBM) has led to an extensive literature on QBM p(i| ], 
particularly in regard to decoherence |^l|,^5| . 

Ford, Lewis, and O'Connell have done systematic and 
comprehensive work on the problem of charge motion in an 
electromagnetic field as a thermal bath in the linear (dipole 
coupling) regime p6|] . They have detailed the conditions 
for causality and a good thermodynamic limit, and have 
further used the QLE paradigm to find stochastic equa- 
tions of motion for non-relativistic charged particles in the 
equilibrium limit. A crucial point of their analysis is that 
the solutions can be (depending on the cutoff of the field 
spectral density) runaway free and causal in the late-time 
limit [^7| . In [£8| , they suggest a form of the equations of 
motion that give fluctuations without dissipation for a free 
electron, but this result arises from the particular choice 
they made of a field cutoff. In contrast, we take the ef- 
fective theory point of view which emphasizes the typical 
insensitivity of low-energy phenomena to unobserved high- 
energy structure. The following two papers in this series 
make clear that a special value for the cutoff is unnecessary 
for consistency of the low-energy behavior^. Further dis- 
tinctions are briefed as follows: First, our general method 
extends to the nonequilibrium regime, allowing considera- 



2 Though there is a maximum value of the cutoff beyond which 
the dynamics become unstable. 
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tion of specific initial states at a finite time in the past. 
Second, we do not make the dipole approximation ab initio 
because we are especially interested in the nonequilibrium 
dynamics of nonlinear systems. When our equations of 
motion are linearized, we pay special attention to the con- 
straints that the original nonlinear theory imposes. Third, 
by adopting the worldline quantization framework, we are 
able to derive equations of motion from fully relativistic 
quantum mechanics. Fourth, we highlight the role of deco- 
herence in the emergence of a stochastic regime. Fifth, we 
show how nonlocality may still be a feature in the stochas- 
tic regime, unless sufficient noise washes out the nonlocal 
correlations between separate particles. 

Our results also represent an extension of Barone and 
Caldeira's analysis of decoherence for an electron in an 
quantum electromagnetic field |50j. Their work employs 
the path integral method except they are limited to non- 
relativistic particles and dipole coupling, and it focuses 
mainly on the reduced density matrix. Rather we give a 
full description of the stochastic dynamics of this nonlocal 
nonlinear theory. An advantage of Barone and Caldeira's 
work is that it is not limited to initially factorized states; 
they use the preparation functional method which allows 
the inclusion of initial particle field correlations. Despite 
this, Romero and Paz have pointed out that the prepara- 
tion function method still suffers from an unphysical depic- 
tion of the initial state |3^]. For this reason, a completely 
satisfactory treatment of particle decoherence in a field- 
environment that accounts for particle-field correlations in 
more realistic (i.e. physically prepared) states remains to 
be given. 

Using the influence functional, Diosi pl| derives a 
Markovian master equation in non-relativistic quantum 
mechanics. In contrast, it is our intent to emphasize the 
non-Markovian and nonequilibrium regimes with special 
attention paid to self-consistency. The work of |32j dif- 
fers from ours in the treatment of the influence functional 
as a functional of particle trajectories in the relativistic 
worldline quantization framework. Ford has considered the 
loss of electron coherence from vacuum fluctuation induced 
noise with the same noise kernel that we employ [ p3[ . How- 
ever, his application concerns the case of fixed or predeter- 
mined trajectories. 

F. Organization, notations and units 

In Section 2, we begin with a review on how a quan- 
tum field can be treated as a bath of harmonic oscillators 
and we connect with the well-studied quantum Brownian 
motion model (QBM) and quantum Langevin equations 
(QLE) 0^g,|||g,|l3|. We write down the form of 
the influence functional (IF) for nonlinear particle-field in- 
teractions, and we further review the influence functional 
formalism in Appendix A, where we discuss the necessary 
assumptions for deriving an evolution propagator for the 
kind of hybrid particle/field model that we employ. Under- 



lying the formulation in Appendix A is the worldline quan- 
tization method for relativistic particles where the particle 
coordinate is the relevant particle degree of freedom. Use 
of this framework, rather than the conventional quantum 
field description for the charged particles, is central to this 
work. However, the detailed development of this approach 
is not needed for the semiclassical and stochastic regimes 
considered in this first series. 

In Sec. 3, we derive the coarse-grained effective action 
(CGEA). The background material in our approach con- 
sists of the Schwinger-Keldysh closed-time-path (CTP) ef- 
fective action applied to an open system, which results in a 
CTP CGEA that is closely related to the influence action in 
the Feynman- Vernon influence functional method. For easy 
access and identification of notations and conventions, we 
review these necessary formalisms in the Appendices, which 
can be read independently. These techniques provide a gen- 
eral and powerful framework for studying nonequilibrium 
quantum processes, especially for non-Markovian dynam- 
ics, which are prevalent when backreaction of the environ- 
ment on the system is fully and correctly accounted for. 
We introduce the reduced-density-matrix for a system af- 
ter assuming a system-environment partition. The reduced 
density matrix evolution operator is found in terms of the 
IF, and the CGEA. In Sec. 4, we define the stochastic ef- 
fective action, and show the close relationship of noise, dis- 
sipation, and decoherence in the stochastic limit. In Sec. 5, 
we discuss stochastic fields. Sec. 6 presents a general (i.e., 
not restricted to a specific model of the particles or field) 
derivation of the stochastic equations in the form of nonlin- 
ear Langevin equations for general nonlinear particle-field 
coupling. We discuss when the traditional linear Langevin 
equations may be recovered. In Sec. 7 we derive a gen- 
eralized fluctuation-dissipation relation for the noise and 
radiation-reaction in these expressions. This provides the 
necessary framework of a stochastic theory approach to in- 
vestigate relativistic particle motion in quantum fields. 

In Paper II |E9[ , we shall apply these results to a specific 
model of relativistic particles in a scalar field. We derive 
the influence functional and stochastic effective action for 
this nonlinear model, and then discuss in detail the semi- 
classical limit. The main result is the derivation of equa- 
tions of motion for the semiclassical limit that, to lowest 
order, are modified (time-dependent) Abraham-Lorentz- 
Dirac (ALD) equations for charged-particle radiation re- 
action. This work demonstrates that the ALD equation is 
a good approximation to the semiclassical limit of scalar 
QED, in the regime where the particles are effectively clas- 
sical. This derivation is more general, however, and reveals 
the time-dependent renormalization and radiation-reaction 
typical of nonequilibrium quantum field theory. These fea- 
tures play a crucial role in establishing the consistency of 
this limit by showing that runaway solutions and acausal 
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effects do not occur when the field is suitably regulatedg. 
We then derive (time-dependent) Abraham-Lorentz-Dirac- 
Langevin (ALDL) equations, which describe the fluctua- 
tions (and radiation reaction) of relativistic particles in the 
stochastic limit. We use these results to explore the ques- 
tion of Brownian motion for a free particle in a quantum 
scalar field. 

In Paper III Q, we extend these results to the electro- 
magnetic field by deriving the influence functional for QED 
in the Lorentz gauge, from which we find the interacting 
photon-particle stochastic action. In the semiclassical limit 
we obtain a modified, time-dependent, ALD solution for 
QED. We then demonstrate a direction dependent Unruh 
effect with the vacuum fluctuations of the field appearing 
thermal (at the Unruh temperature) for charged particles 
in a constant electric field. This general approach clearly 
illustrates the distinct natures of ALD radiation reaction, 
the associated ordinary radiation into infinity, and the Un- 
ruh effect, while at the same time deriving all these effects 
in a unified and self-consistent manner. As an interest- 
ing application, the Langevin equation may be viewed as a 
model for a stochastic particle event horizon, with implica- 
tions for a nonequilibrium treatment of black hole Hawking 
radiation and horizon fluctuations. 

Our second series of papers (IV and V |3(|) will be on 
relativistic quantum particle-field interaction. It is here 
that the pre-cursors of the quantized worldline framework 
introduced in Appendix A are fully developed. In Paper 
IV, we 'first' quantize the free relativistic particle using 
the path integral method. The worldline path integral rep- 
resentation is used to construct the "in-out" and "in-in" 
relativistic particle generating functional for worldline co- 
ordinate correlation functions. In Paper V, we use the 
influence functional to construct an interacting quantum 
theory of both relativistic particles and fields. This in- 
volves the introduction of a nonlocal worldline kernel that 
provides a route for exploring the role of correlation, dis- 
sipation, and nonequilibrium open-system phenomena in 
particle creation and other relativistic processes. This work 
is also a model for how open-system methodology may be 
introduced into String theory since the worldline path in- 
tegral is the point particle limit of the String case. This 
formalism provides a powerful approach to a new set of 
quantum-relativistic-statistical processes where one may si- 
multaneously take advantage of the influence functional 
and first-quantization (worldline) techniques. 

We shall usually set c = 1, but keep H explicit as a 
marker for quantum effects. The metric tensor is 
diag(+l, — 1, — 1, — 1). Greek letters will (usually) denote 
spacetime indices, and Latin letters a, 6, c will indicate 



3 This demonstration of causality addresses the early-time, 
nonequilirium setting rather than the late-time equilibrium limit 
that is analyzed in At late times, the conclusions in ^] 

also apply to the linearized version of our results. 



the CTP time-branch (see Appendix C). The mixed func- 
tion/functional notation f[z;x) will be used when / is a 
functional of z, but a function of x. The Einstein summa- 
tion convention is employed except when otherwise noted. 
Particle degrees of freedom will be collectively denoted 
z = {z£ (t„)} with r„ as their worldline parameters. The 
path integral measures are denoted Dz and Dip, and are 
defined in Appendix A. 

II. THE EFFECTIVE ACTION, EFFECTIVE FIELD 
THEORY, AND PARTICLE REDUCED 
DENSITY-MATRIX 

Particles moving in quantum fields may be cast in the 
form of nonlinear quantum Brownian motion. In this pa- 
per, a set of N particles (with action S z ) constitutes the 
system and a scalar field ip (x) (with the action S v ) con- 
stitutes the environment. The particles and field interact 
nonlinearly as determined by an interaction action Si n t- 
The n th particle degree of freedom is its spacetime trajec- 
tories z% (t„) parametrized by r„, which is not necessarily 
the particle's proper time. In this first series, our focus is 
on the stochastic dynamics of the particle worldlines; for 
simplicity, we neglect additional particle degrees of freedom 
such as spin. We also work in the regime where particles 
have approximately well-defined (emergent) trajectories, 
though with statistical fluctuations arising from the par- 
ticles' contact with the quantum field. Consequently, the 
forces that arise from spin-statistics (e.g. Fermi-exclusion 
or Bose-condensation) are negligible. We discuss how these 
quantum effects are suppressed by decoherence (the same 
decoherence that gives the particles stochastic trajectories) 
in our second series. 

In this relativistic treatment, both the particle's space 
and time coordinates will emerge as stochastic processes, 
with correlation functions of the form (z^ (n) ...z v N (tn)) s . 
The notation () s denotes the stochastic average with re- 
spect to a derived probability distribution. The statistics 
of these processes are found from the underlying quan- 
tum statistics of the quantum field acting as an environ- 
ment. The semiclassical regime is characterized by the par- 
ticle trajectories being well-defined classical variables. The 
stochastic regime is intermediate between the semiclassi- 
cal and quantum regimes, where the particle worldlines are 
classical stochastic processes. 

The quantum open-system of particle degrees of freedom 
is described by the reduced density-matrix 

Pr (t) = Tr v p v (t) , (2.1) 

where pu (t) is the unitarily evolving state of the universe 
(i.e. the particles plus field) at time t. We assume that the 
initial state (at time ti) of the universe can be represented 
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in the factorized (tensor-product) form 

pu[ti) = Pz{U)®P v [ti) (2.2) 
dipidifiidztdz'tp^ (tpi, 

X p z (Zi,z'i) \<fi,Zi)((Pi,zl\, 

where p z and p v have been expanded in terms of basis 
states | ipi, Zi) — \ipi}® |zj, ti). The initial particles states are 
the Lorentz invariant relativistic configuration-space state 
defined as 



I*) 



d 4 P 



(2.3) 



(2-r 



p — m 



>°) 



\P), 



where we have assumed that the initial particle state is 
positive frequency (this assumption may be relaxed). 

The particle's reduced density-matrix at later times is 
given by (see Appendix A) 

p r (zf,Zf)= J dzidz' { J r (zf , z' f ; Zi, Zi) (2.4) 
X Pz {zi,z[) . 

The variable z stands for the entire collection of particle 
coordinates {zfi} . The open-system evolution operator J r 
is given by 



J r [Zf, Zf \ Zi, zl) 



DzDz'e 



>MsAz]-S*[z']) 



x F [z, z] 



(2.5) 



F [z, z'\ is the influence functional introduced by Feynman 
and Vernon Q. The measure Dz depends on the details 
of how the particle action is defined and quantized, and it 
may implicitly include additional gauge degrees of freedom 
(such as a lapse variable N). In the semiclassical/stochastic 
treatment presented in this paper, we will in fact assume 
that the particle degrees of freedom have been gauge-fixed 
so that the parameters r are the proper-times. The detailed 
treatment of the influence functional for relativistic particle 
worldlines is found in our second series p(| , though a brief 
discussion may be found in Appendix A and in paper III 

The influence functional is given by 

r r<Pf=<p' f 
F[z,z']= dipfdtfidip', / DipDip' (2.6) 

X Pip (<Pi,<Pi',ti) 

= exp j^SV-F > 

where Sjf [z, z'\ is called the influence action. Notice that 
the definition of the influence functional itself involves only 



path integrals over the fields <p, and hence, does not require 
doing particle worldline path integrals. Indeed, the results 
of this paper follow from effectively treating the worldlines 
as classical variables coupled to a quantum field, which is 
the usual semiclassical paradigm. In the next section we 
indicate the connection between the worldlines as classi- 
cal variables and the quantum equations of motion for the 
quantum-average trajectories. 

Defining the coarse-grained effective action 



ScGEA = S z [z] ~ S z [z'] + S IF [Z, Z 1 ] . 

J r may, alternatively, be expressed as 



[z f ,z f ; 



DzDz'e^ ScGEA \- z 



(2.7) 



(2. 



The evolution kernel J r gives all information about the 
quantum open-system dynamics of the particles. 

If the particle-field interaction is characterized by a cou- 
pling constant e and is linear in the field variables, the 
influenc e ac tion is of order e 2 . One approximation scheme 
for Eq. (2.8) is to expand the influence functional as 



(2.9) 



giving a perturbation theory based on the order of a cou- 
pling constant. J r then describes free quantized-particle 
evolution punctuated by discrete interactions with the field. 
This is just ordinary perturbative quantum field theory in 
the "in-in" formulation. Unlike "in-out" field theory, the 
lowest order (non-trivial) terms are of order e 2 , which is 
natural since J r gives the evolution of the density matrix, 
not the wavefunction. But , p articularly when there is a 
background field (p = (<p), (2J3) does not give an accurate 



depiction of the classical limit of particle trajectories, say, 
as given by the mean- values (z n ) = z, where z n is the quan- 
tum operator for the particle's coordinates. The infinite 



set of background field terms in the expansion from (2.9) 
may be summed, giving a far more accurate approximation 
to the particle evolution. The resulting re-organization of 
Feynman diagrams is automatic ally obtained by perform- 
ing a saddle point expansion of (2 



It is convenient to define sum and difference variables 



g~ = (z- z') 
z+ = (z + z')/2. 



(2-10) 



Evaluating the evolution kernel in Eq. (2.8) for linear 
systems is a standard step in the derivation of the mas- 
ter equation using path integral methods. One begins by 
reparametrizing the worldlines defining the fluctuation and 
classical worldlines, z n and z n , respectively, by 

^ = ^-^> (2-11) 

The extremal solutions to the effective action, denoted z n , 
are defined to be the classical solutions to the real part of 
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ScGEA [z d 



= o, 



satisfying the boundary conditions 



(2.12) 



(2.13) 



This definition is made because the imaginary part of the 
effective action does not modify the stationary phase solu- 
tion. 

The evolution operator for the open system is then given 

by 



J r (zf,zf) = e^ s cGE A [^\ J J D z± (2.14) 

, 5 2 S r 



x e 



x (l + 0(z 3 )), 

where we have factored out the "zero- loop" term in front, 
and expanded the exponential of the non-quadratic terms 
in the effective action. S GGEA is the imaginary part of the 
effective action. 

The higher-order quantum corrections depend on the 
solutions z„ unless the effective action is quadratic; and 
hence, the model is purely linear. For the linear case, J r 
may be found exactly since all the path integrals are Gaus- 
sian)^]. It is given by 

Jr (zf\z?) =F{t f ,ti)x eiSoasAis*] (2.15) 

where F may be found by the normalization condition 



dzfdz' f 6(z f - z' f ) J r (z f ,z' f ; zf |^ =i , 
S(zi-zi). (2.16) 



From J r it is straightforward to find the master equation 
from the expression 

^-p r (zf) = Jdz+dzrJL Jr (zf;zf) p r {zf) . 

Because the one-loop (fluctuation) term is independent of 
z for linear theories, the resulting noise and dissipation 
are independent of the system history. It is only for lin- 
ear theories that the noise induced by the environment is 
completely extrinsic in this sense. 

For non-linear theories, the one-loop (and higher-order) 
corrections will depend on z^ 1 and hence will modify the 



4 Actually, the integrals are nearly, but not quite, Gaussian 
because the limits on the time-integrals is not [— oo,+oo] but 
rather [ti, +oo] , where ti is the initial time. 



master equation, giving system-history dependent (col- 
ored) noise and dissipation. But, rather than using the 
master equation as our primary tool, we instead follow the 
approach in |],|37| for deriving quantum Langevin equa- 
tions EH], to find a stochastic description of the particle 
dynamics from the influence functional. 



III. INFLUENCE FUNCTIONAL FOR QUANTUM 
SCALAR FIELD ENVIRONMENT 

We now derive the influence functional assuming a sys- 
tem of spinless particles locally coupled to a scalar field. 
The free scalar field action is 



Sb[<p] 



dt dx 



1 



(3.1) 



We assume that the interaction term is linear in the field 
variables <p : 



Sint = dtdx j[z n (t„) ; x,t) ip(t, x), 



(3.2) 



but j[z;x,t) is an arbitrary (nonlinear) functional of the 
particle coordinates z n . The construction of a perturba- 
tion theory to derive the influence functional for nonlinear 
environments (e.g. y 4 ) has been developed by Hu, Paz and 
Zhang in [ pT| . We consider couplings that are nonlinear in 
the system variables, as is the case for QED. For a sys- 
tem interacting with a quantum field, one can follow the 
method introduced by Hu and Matacz |34|] for deriving the 
influence functional in terms of the amplitude functions of 
the parametric oscillators which are the normal modes of 
the field. We expand the field in terms of its real-valued, 
normal modes 



WAfc(x) 



as 



sink ■ x ; A = 1 
cosk • x; A = 2 



tp(x) - (2/L) 3/2 £ u kA (x) 

kA 

= (2/L) 3/2 J>aW Ua(x), 



-\/k2 



(3.3) 



(3.4) 



where a = (k, A) . We assume that the field is in a box 
of size L with periodic boundary conditions so that the 
mode wave vectors are given by k = 27rn/L for positive 
integers n^. We may take the limit L — > oo to recover the 
continuum limit, though in doing so we must be mindful of 
both Lorentz invariance and infrared divergences. 

With this normal-mode decomposition, the environment 
becomes a collection of real harmonic oscillators ip a (t) that 
are linearly coupled to the system with the quantum Brow- 
nian motion action 

SB + S int = J2 f dt [{<p 2 a /2 + wl<f%/2 + <p a f a (t)}, 

a J 

(3.5) 
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where 



f a (t) = (2/L) 3/2 / dx i[«(t);x,t)u a (x). 



(3.6) 



The functional integrals (in the path integral representa- 
tion of the influence functional Eq. ( |2.6| )) then become 
integrals over the oscillator's coordinates. The Jacobian 
for the change of measure J ±I X t dcp(t, x) — > J II Q d</? a (i) 
is 1. We may therefore use the well-known results for the 
influence functional |H]|34j ■ Initially correlated states may 
be treated by the preparation function method [p2f , but the 
qualitative results are largely the same as those for initially 
uncorrelated states after fast transients have subsided [[5lJ . 

We note that it is a consequence of the fact that the in- 
fluence functional is derived in the functional Schrodinger 
picture that there is an associated choice of reference frame 
and fiducial time. Hence, even in the L — > oo limit, the 
influence functional is not manifestly covariant, but in- 
stead depends explicitly and implicitly on the initial time 
at which the particle and field states are defined. This ini- 
tial time is also identified with the spacelike hypersurface 
on which the initial state takes a factorized form, which 
is likewise not a covariant construction (See appendix C 
for more discussion on this point). Despite these facts, we 
see in Paper II that the equations of motion for the par- 
ticles quickly 'forget' the initial time U, and therefore be- 
come well-approximated by manifestly covariant equations 
of motion at later times. 

The influence functional factorizes as a product of terms 
F a for each field mode a : 



=exp{~S IF [f-,f+]j=Il a F a [ft] (3-7) 

= cx p{-^E / dt f dt'{2f~(t)p a (t,t>)f+(t>) 

-\f-{tW{t,t')f-{t')]}, (3.8) 
We have defined sum and difference variables 



#(*) = [/<*(*) + £(*)] /2 

f«(t) = [f a (t)-fL(t)}, 



(3.9) 



in terms of which Sif takes a particularly simple form. 
The influence kernel 

Ca(t,t') = v a {t,t')+ip a [t,t') (3.10) 

for a field initially in a thermal state 

p ip =exp{-0H v \ (3.11) 

(at inverse temperature (3 = l/ksT) is given by 

£„ = {coth (h(3u) n /2) cosuj n (t — t') — ismu) n (t — <')}. 

(3.12) 



Substitution of ( p.6[ ) into ( |3.8| ) allows us to write the influ- 
ence functional as 



F\j(z),j(z')} 

= exp"^-^— J dx J dx' (j[z;x)~j[z';x))x 
(G + (x,x')j[z;x') + G + (x,x')*j[z';x'))\ 



= exp \ -Si F [j,j'} 



(3.13) 



(3.14) 



where G + (x,x') is the positive frequency Wightman func- 
tions given by the mode expansion 



G + (x,x') = 




: °tt> U a (x) U Q (x')Ck(t, (3.15) 



cos k(x-x')Ck(*,*')- ( 3 - 16 ) 



In going from ( |3 .15 ) to (3.16) we have 
used ^2 X «kx(x) itkA(x') = cos k(x — x'). We abbreviate 
9 (t, t') = 6 tt i . The factor of 1/2 comes from double count- 
ing the equivalent modes k and — k. In the limit L — > oo, 
we may go from discrete to continuous modes by replacing 

Sk>n -» (L 3 /2) 1/2 f °° dk. The Green's function G+ (x, x') 



is given by the free field correlation function 

G + (x,x')=Tr(ip(x)ip(x')p v (ti)) 



(3.17) 



where p v (ti) is the initial state of the field. The real part 
of G + is the field commutator (also called the Schwinger, 
or causal, Green's function G c ): 



G ( 



ReG+(x,x') = ([<p(x),<p(x')]), 



(3.18) 



The imaginary part of G + is the field anticommutator (also 
called the Hadamard Green's function): 



G 



H 



JmG + (x,x') = ({tp(x),cp(x')}) 



(3.19) 



The commutator is a quantum state independent function; 
it determines the causal radiation fields in both the quan- 
tum and classical theory, and is responsible for dissipa- 
tion in the equations of motion. G H is quantum state de- 
pendent; it determines the form of quantum correlations, 
and is responsible for noise in the stochastic limit. In- 
deed, because the field-part of the action is quadratic, the 
Hadamard Green's function contains the full information 
about the field statistics; though (as we shall see) the field- 
induced noise is determined by the quantum field statistics 
plus the particle kinematics. 

Using the symmetry properties of the Green's functions 
and the form of the influence functional, only the retarded 
part of the commutator appears in the equations of motion, 
where 

G R (x,x') = 9(t,t')ReG+{x,x'). 
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Defining sum and difference variables 

j+ = \j[z,x)+j[z',x)]/2 
j~ = \j[z,x)-j[z',x)] , 

the influence functional takes on the form 



(3.20) 



^[J't']= ex p{x J dx J dx'[23~(x)G R (x,x')3 + (x>) 

-ij-(x)G H (x',x')j-(x')}}. (3.21) 

When the function j(z,x) factors as j(z, x) = 
z n h n {x), we recover the influence functional for lin- 
ear quantum Brownian motion of n particles follow- 
ing arbitrary yet prescribed (i.e. not dynamically nor 
sclf-consistently determined) trajectories. The functions 
h n (t,x) provide effective time-dependent coupling con- 
stants |54j , the influence functional for this case is reviewed 
in Appendix B. The closely related CTP coarse-grained ef- 
fective action, and its adaptation to the particle-field mod- 
els we consider, is reviewed in Appendix C. Further use of 
the CTP formalism, together with the nonlinear particle- 
field influence functional, is made in the second series, 
where the quantum regime of relativistic particles is de- 
veloped. 

A. The stochastic generating functional 

Rather than developing the the ory at the lev el of the 
master equation, we return to Eq. (2.4) and ( |2.8| ) but now 
add the new source terms 



dr(j,(r) Z " (r)-J» (r) ) 



dT{J+z" 



(3.22) 



to the coarse-grained effective action (CGEA), set Zf 
and integrate over Zf (with the restriction 
Appendices). We find 



> U 



see 



Z[J±] = J dZfPr(Zf,Zf)j± 



f f z f~ z f 
= / dzfdzidZj / DzDz' 

x e i(S C GEA[z,z'] + f (j+ g~ +.J- z + )dr) 

xp(zi,z'i), (3.23) 

where Z [J^] is the precisely the generating functional 
shown in Appendix C on closed-time-path methods. Z [j^j 
may be used to derive the particle correlation functions 



(z (t)), (z (t) z (t')), ... via , e.g., (C31), or as the starting 



point for deriv ing the 1PI coarse-grained effective action 
via, e.g., ( |C32| ). In this section, to keep expressions as sim- 
ple as possible, we will suppress most indices (e.g. Lorentz 



indices, particle labels) and return to more explicit nota- 
tions later. We also won't worry about any overall constant 
normalization factors since these don't effect the resulting 
equations of motion. 

The crucial observation for translating a stochastic de- 
scription hinges on the fact that the noise kernel (appearing 
in ImSip) is a real, symmetric kernel with positive eigen- 
values. Therefore, 



g ~k ScGEA 



= P -iiSj F = P 4EJ T dtdt' f-(t)u a (t,t')f-(t') 



= \F[z,z']\<0. (3.24) 

The inequality holds for any pair (z, z') . Recalling that 
/ _ = f[z] — f [z 1 ] and z~ = z — z' , we may expand the 



argument of the exponential (3.24) in powers of z giving 



ex P {-i^y dtdt'[(Sf a (t)/dz^(r))^ a (t,t') 



x (6f a (t>)/6z»(T>)) 
+ (.-■)}. 



z-=0 



z*- (r) z»- (r') drdr' 



(3.25) 



The influence functional is exponentially small for "large" 
deviations of the histories z~ (r) . Of course, what counts 
as large depends on the noise kernel and coupling repre- 
sented by / it) . (This suppression of the magnitude of the 
influence functional is indicative of decoherence. Off diag- 
onal terms in the density matrix p ( z~^ , zj J will have large 



z and hence tend to be suppressed by (3.25).) 



On the assumption that decoherence is strong enough 
(whic h nee ds to be verified for particular cases and mod- 
els), ( |3.25 ) justifies the expansion of the effective action in 
powers of z~ . The real part of the CGEA (including source 
terms) is 



dT 



5z^- (r) 



ScGEA \ z i z '\ 



+ J+(t))z»- (t) 



+ f drJ-(r)z^+ (t)+0(z- 3 )]. 



(3.26) 



Together with ( |3.25| ) , the generating functional is (neglect- 
ing O (z~) 3 terms) 



Z[J±] 



(3.27) 



dzf dz 



Dz + Dz~ 



x exp 



dr 



SS. 



CGEA 



z~=0 



I, ^ J \ <):■' (V) 

( T ) +i Yl J dt ' dTdT ' 



z~=0 



\5z»+ {t)J " v ' ' \SzK+ [t') i 
z»- (r) z»~ (r') + / drJ~ (r) z^ (r) }p (zf) 
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In general, the initial state of the particle p (zf) is arbi- 
trary, but we are interested in finding a description of the 
particle's stochastic trajectory when the particle position is 
well localized. We therefore take the initial particle density 
matrix to be 



Pz [ti) = \Zi,ti){Zi,ti\. 



(3.28) 



The generating functional takes a simpler form in this case; 
we find 



Z[J±] 



(3.29) 



= / dzfdz 



x exp 



Dz + Dz~ 



CGEA 



x z r 



(r)+iJ2 J dt'drdr' 



8z^- (r) 



z»-(t)z»- (r')+ / drJ-(T)z^{r)}. 



In our second series of papers we discuss more general ini- 
tial states. 

The manipulations above are somewhat formal in that 
we have not explicitly discussed how to perform the z ± 
path integrals (discussed in our second series). These de- 
tails are important when considering higher-order quantum 
corrections, but for the semiclassical/stochastic limit it is 
enough that we can do standard Gaussian path integrals 
in z , which we can. In the Gaussian approximation the 
only important modification comes from the restriction of 
the integration range placed on the particle time variable: 
ti < z <tf (see Appendix A). In the limit that initial and 
final states are defined at ti — — oo and tf = +00 (as is the 
case for asymptotic scattering processes) the Dz° integra- 
tions are the usual Gaussian- type. More generally, when 
considering processes not too close (in time) to the initial 
time hypersurfaces X (ti) , the restriction on the z° inte- 
gration range has negligible effect on Z [J*] because the 
integration range of z° may be extended to ±00 with negli- 
gible error. But, for very short times after ti, the generating 
functional is modified by the z° boundary conditions. This 
has the result of modifying both the noise and dissipation 
in the Langevin equation derived below at very short times 
after ti in such a way as to preserved consistency with the 
given initial state. 

These modifications are in the nature of what O'Connor 
Stephens and Hu refer to as finite-size effects arising from 
the existence of boundaries or other topological (or more 
generally even dynamically-generated) constraints in fields 
or spacetimes |38| . More familiar are finite size effects com- 
ing from spatial boundaries, implying a spatial restriction 
on the paths in the sum over histories. Periodicity in (imag- 
inary) time depicting thermal behavior may also be viewed 



as a kind of finite-size effect . One advantage of the space- 
time formulation (in terms of worldline path integrals) is 
that it allows a natural description of finite-size effects aris- 
ing from both spatial and temporal (i.e. spacetime) bound- 
aries or restrictions. Hawking and Hartle's path integral 
calculation of black hole radiance may be viewed as a fi- 
nite size effect where the singularity and associated black 
hole geometry restrict the paths leading to thermal parti- 
cle creation at the Hawking temperature [^9| . Similarly, 
Duru and Unal |l0| use path integrals to calculate particle 
production in a cosmological spacetime where the initial 
singularity places a temporal restriction on the particle 
paths similar to the restrictions giving in this paper. In 
p0| , the restriction on worldlines (z° (r) > ti) is required 
by the initial cosmological singularity at ti, in our work 
the restriction is part of the initial condition: we take the 
particle state to be known (e.g. fixed) at, and before, ti. 
Consistency then requires the restriction z° (r) > ti. 
Defining the variable 



^ it) 



6z»- (t) 



S(JGEA [ z± ] 



+ J+ (t) , (3.30) 



we may do the (approximate) Gaussian Dz path integral 
in Q3.27|). The result is 



Z [J*] = J Dz+ exp { - i J J drdr' 
x, 7m (t) C-^t.tO/M 
+ y drJ-(r)z^(r)} (3.31) 

The kernel C _1 is defined by 



5 (r - r') = / d T " dt!' dt! C~ x (r, t") x 



(3.32) 



5z»+(t)J \6z+(t') 

= \ dT"C- X {T,T")C(T" ,T') 



We denote the kernel C (the inverse of C- 1 ) by to 
indicate its association with rj (r) , and to dis tingu ish it 
from other kernels that arise later. Next, Eq. j3~30| ) may 
be solved for functional of 77 (r) and J + (r) (i.e., 

z + = z + [rj, J + ])0- We now make the change of integration 
variables Dz + — > Dr\ using the fact that the Jacobian is 
equal to 1. Then the term 



P [77] - N a e 



J dTdT , Vlt (r)C- 1 (r,r') V "(Tt') 



(3.33) 



5 In the more common linear case, the noise is independent of 
the system histories (z + ). This is a crucial difference between 
linear and nonlinear theories. 



12 



may be interpr eted as a probability distribution for a noise 
rj (t) , turning (3.30) into a Langevin equation. Nq is some 
appropriate normali zation constant. 

In the expression (3.33) above, the kernel C^ 1 is a func- 
tional of z + . This fact makes the noise system- history de- 
pendent (non-Markovian). With the change in path inte- 
gration variables, we find 



Z J d 



D V P [ T j\etf* rJ »W*S +( - T \ 



(3.34) 



where we note that 77 is a functional of J + through Eq. 
(3.30). We have added the subscript 77 to z + to indicate 
that the change of variables may be inverted to give 
as a functional of 77 via the Langevin equation (3.30). For 
simplicity, we set J + = 0, and write z + — > z and ,/ + — > J. 
The generating functional 



Z[J] = J DrjP [77] / H^KW 



(3.35) 



may then be used to find the trajectory correlation func- 
tions. The mean-trajectory is 



ft 6Z[J] 



i <5J M (r) 



= / Dr]P [rj\ [<(r)] =<^(r 1 )) 7) 

,7=0 J to 

= z» (r) . (3.36) 



The n-point correlators for the particle are given by 
ft™ S n Z [J] 



l n 5J li (T 1 )...5J v (T n ) 



(3.37) 



j=o 



Dr,P[n} [z* (n)...< (r„)] 



= <*" (t x )...z v (r n )) v . 
Brackets with a subscript 77 denote the stochastic average 



(y(ri)) ( 



Dr)P[r)]y (77) . 



(3.38) 



Since C 1 in ( |3.33 ) is still a functional of rj (through Eq. 
(3.30)), this is not a Gaussian probability distribution, but 



is far more complicated involving the full effect of backre- 
action between the system and field. 

To mak e pro gress, we linearize the nonlinear Langevin 
equation ( 3.30| ) by expanding the trajectory around the 
solution, z (t) , to 



5z^- (t) 
5 



Sz»- (r) 



ScGEA [ z± ] 



ScGEA [z ] 



(3.39) 



=0,z+=z 



:-=0,z + = 



= 0. 



This solution we call the semiclassical trajectory; observe 
that the imaginary part of Scgea plays no role at this 
level. Defining 



(3.40) 



wc find 



?7m ( t ) 



SScGEA [z^] 



5z^- (r) 



z —Q.z + —z 
SScGEA [Z^] 

5z v + (t') V Szv- (r) 



+ / dr'z" (/) (3.41) 



+ 0(z 



z2\ 



Since the first term on the LHS of (3.41) vanishes by defi- 
nition of z, we find the linearized Langevin equations 



Vli (r) = / dr'r (r') 



/ 8S, 



CGEA [Z" 



Sz v + (t') V 5z^- (t) 



z~=0 

(3.42) 



Expanding the probability distribution around z gives 

P[ V ] = N e-&ff drdr'^(r)C^{ry)^(r')+0{^)^ (3 ^ 

since 

C" 1 (r, r') = CJ 1 (r, r') + (SC^/Sz) S + O (z 2 ) , 

and z is O (77) . At this point, one can make a Gaussian 
approximation to the noise, which involves evaluating C _1 
with the (self-consistently determined) semiclassical solu- 
tions z. Later we shall show how to treat non-Gaussian 
noise via a cumulant expansion. 



Under a Gaussian approximation, it is clear from ( 3.43 ) 
that the mean noise (77)^ = 0. The noise correlation func- 
tion becomes 



Mr) Vv [t')) = 2 9ijiv C^{t : t i ) 



B. The stochastic effective action 



(3.44) 



From the preceding section we are able to justify a sim- 
pler method of deriving the stochastic equations of mo- 
tion from the stochastic effective action. We also see more 
clearly the underlying role of the quantum field in produc- 
ing noise in this approach. Returning to the expression for 
the influence functional in ( p. 24 ), we use a standard func- 
tional Gaussian identity to write the imaginary part of the 
influence action as 

\F[z,z']\ = N J D£exp{ - J £ (3.45) 

Q 

x ( J dtdt'Ut^^t^ut') 
+ l -J dtut)f*(t))} 

= J £>£P[£]e X p J dtUt)fa (*)| 

= |expj^ J #£ a (*)/a(*)j^ , 
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where £ Q (i) are stochastic variables, and 

P[£] = AT exp j--^ J dtdt'Ut)v«Ht,t')Ut')^ 

(3.46) 

is a normalized probability distribution on the space of 
functions {£, a (t)} . The kernel v~ x is defined by 



(3.47) 



dt'v- x {t,t l )v ot {t l ,t") = 5{t-t"). 

This allows us to write 

exp \i (ScGEA ^CGEA )/n] 



exp( i:Sg GEA ) \F[z,z'} 



(3.48) 



as 



e 1cSoG*A[*,*'] = J D ZP[Z]expl_{s A [z] - S A [z'\ 

+E / */«(*)(&(*) 

+ 2 f dt'ii a (t,t')f +(?))} 
= J D£P[£]e^ s *[ z ' z 'L (3.49) 
This expression defines the stochastic effective action 
S t [z, z'\ = S A [z] - S A [z'\ + W # /a (*)(&(*) 

+ 2 / dt'^ a (t,t')f+(t')). (3.50) 



The variables £ a (t) are stochastic forces (noises) coupled 
to the currents f a (t). The mean force vanishes since 

(&(*)>« (3-51) 
= iVo / D^eaWe { "^ /dtdt ' e " (t) ^ 1(i ' t,)ea(t ' )l} =0. 



The £ noise correlator is 

({&(*), &(f)}> € = J D^ a Ut)Ut')P[U 

= hv a (t,t'). (3.52) 

Notice that in these expressions no Gaussian approxi- 
mation to the noise is made. The noise variables £ Q (t) , 
unlike the earlier noises T] (r) , are already exactly Gaus- 
sian because their statistics are precisely those of the (as- 
sumed) Gaussian quantum field state. On the other hand, 
the stochastic noises £ a (t) do not couple directly to the par- 
ticle coordinates z (r) , as was the case for the r\ (r) , but 



instead couple to the nonlinear function f a (t) . We shall 
explore this connection in more detail below. 

We may now find the nonlinear Langevin equations of the 
previous section directly from S% as the stochastic solutions 
to 



Sz^- 



0. 



(3.53) 



z~=0 



with the noise characteristics given by (|i5|). As 
in ( 3. 39] ), the semiclassical solutions are found from 
{SScgea/Sz'- 1 ^) z ^ =q — 0. The solutions, being derived 
within a self-consistent initial value ( "in- in" ) formalism, 
are real and causal [^lj. The stochastic regime manifests 
at a deeper level than the semiclassical limit. In our par- 
ticular case, the stochastic regime identifies the quantum 
field statistical correlations with fluctuations induced in the 
particles' trajectories. These fluctuations represent physi- 
cal noise when the particle histories decohere^]. 

C. The cumulant expansion 

We shall find later that higher order effective noise cu- 
mulants arise due to the nonlinearity of the particle-field 
coupling- we saw the first sign of this when we derived the 
stochastic generating functional. To treat this more gen- 
eral case it is convenient to define a generating functional 
for normalized correlation functions 

W[f±] =InF[/±] =(i/K)S IF . 

The cumulant expansion method was used by Hu and Mat- 
acz J37J in their derivation of the Einstein-Langevin equa- 
tion in stochastic semiclassical gravity Ta ylor expanding 
W [f^] in powers of /~ (t) (defined in (3.6)), we find 



W[f+ , /"] = ~ jf dhf- ( tl )c[ f J (h ; /+] 

-^J^dhJJdhf-^f-ih) 
-.(/)/ 



(3.54) 



xCg>(t 1 ,t 2 ;f + } + 

xCW(t 1 ,...,t n ;f+] + 



dt 1 ---dt n f-(t 1 )---f-{t n ) 



6 Whenever we speak of particle histories in the context of the 
semiclassical or stochastic regime, it will be implicitly assumed 
that some additional smearing of the particle histories has been 
done to produce truly decoherent histories (see Sec. 3. 5). When 
the histories are so coarse-grained that they effectively acquire 
substantial inertia (usually the case in the macroscopic limit), 
the efffect of this noise will be comparatively small ]13|. 
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For simplicity, we will drop the mode label a in the re- 
mainder of this Section. The full stochastic action is given 
by the addition of contributions from each mode (since the 
influence functional factors into a product of terms for each 
mode). The real cumulants Cn are given by 



t ■ f H 



5 n W[f+,f-} 



S.f-(h)---Sf-(t n ) 



-=o 
(3.55) 



The C n J > are functionals of f + (t). They are, therefore, de- 
pendent on the system history. The cumulants C n are of 
order e™, with n — 2 being the Gaussian noise contribu- 
tion. The absolute value of the influence functional may 
be written as 



|F[/ + ,r]| = ^exp|| J dtf-(t)E(t)\ 



(3.56) 



#SP[S,/+]exp - / dtf-(t)B(t)}, 



|F[/ ± ]| is therefore the characteristic functional of the 
stochastic process H. The normalized probability func- 
tional P[3, / + ] may be found from the influence functional 
by inverting the functional Fourier transform ( 3.56] )pj The 
probability distribution P[S, / + ] depends on the system 
history / + for nonlinear models, P [£] is independent of 
the / + for linear models. Making a Gaussian noise approx- 
imation by neglecting all cumulants beyond second order 
(i.e. by working to order e 2 in the coupling constant), the 
probability distribution is 



P[S,/+] (3.57) 
= P exp j-^- J dh J (ft 2 3(ti)Ci- 1 (*i>*2;/ + ]3(ta)}, 

and the correlator of the noise S is 

(Z(t 1 )Z(t 2 )) s = hC 2 (t 1 ,t 2 ;f+}. 

In this approximation the stochastic influence action may 
be written 

5^(^,5]= /"dt{/-(t)(Ci(i;/+]+S(t))}. (3.58) 



D. Stochastic fields 

While our primary interest is the dynamics of the parti- 
cle subsystem, it is helpful to see how the quantum field is 



equivalent to a stochastic field in terms of its effects on the 
particle trajectories for the model we are considering, and 
when the stochastic regime is a good approximation. The 
influence functional is derived by integrating out the field 
tp to achieve a description of the effective system dynam- 
ics in terms of particle variables alone together with local 
noise found from the stochastic functions £ a (t). We may 
reintroduce fields by defining the stochastic (field) variables 



% (x)E(2/L) 3/2 ^( a (tK(x), 



(3.59) 



where u a (x) are the normal modes (3.3). The subscript £ 
indicates the functional dependence of x on the stochas- 
tic functions £o,(i). The probability distribution for the 
stochastic fields P[x] is determined by the normalization 
condition 



/fi(x) «/n(0 
where O(x) = { Xs \ £ G 0(0} • Then 

-i -l 



DSP[£], 



(3.60) 



P[xt\ 



But since the change of variables from £ to \ is orthogonal, 
the Jacobian is 1. Therefore, 



P[x] = ex p{ - ^ E J dxdx ' (x(x) Mq(x) 
x v- l {t,t')u a {x!)x{x'))} 



exp { -— / dxdx' x(x) T(x, x') x(x') } , (3.61) 



where 



T(x,x') = ^u a {yi)v a 1 (t,t')u a (-x!). 



(3.62) 



For spacelike separated points (x — x') , T(x,x') does not 
vanish, which reflects the nonlocal character of the quan- 
tum field correlations, and hence the nonlocal correlations 
that will be present between particles that are (even) space- 
like separated. 

The stochastic field x( x ) has vanishing mean, and auto- 
correlation function given by 

(X(x)x(x')) s = K {{${*), = hG H (x, x'). (3.63) 

Thus, x( x ) encodes the same quantum statistical informa- 
tion as the field anticommutator. 



7 In our particular model, the influence functional is Gaussian 
in / , and therefore this inversion is straightforward. 



IV. THE STOCHASTIC EQUATIONS OF MOTION 
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A. General Form in Terms of Stochastic Field 

From the stochastic effective action it is straightforward 
to derive Langevin equations of motion for particle world- 
lines. We emphasize that the exact stochastic properties 
of the noise are known, being derived from the quantum 
statistics of the field- environment. We begin by using the 
definition of \ ( 3.59 ) and f~{t) (3.6) to write the noise 
term in S% as 



E f dt ut)f«{t) 

= I dt f &{j(x,«(t))-j(x ) «'(t))}5]j a (t)tt a (x)) 
dxj~(x)x(x). (4.1) 



solution z!t(r) of (4.4); that is, the system equations of mo- 



Together with the dissipation term from ( 3.21 ), we use (4.1) 
to define 

5 x [z±] = Sa^] + J dxj-(x)( X (x) (4.2) 
+ 2 / dx'G R (x 1 x')j+(x')). 



Thus, the influence of the environment is described in terms 
of a current j~{x) self-interacting through the causal re- 
tarded Green's function G R (x,x'), and with j~ coupled 
multiplicatively to the stochastic field x( x )- 

Using the generating functional W, we find the cumu- 
lants of the noise 



(x\ , . . .X n \ j 



S n W[j+,f 



8j (xi)...Sj (x n ) 



j-=0 



(4.3) 



For the field state initially Gaussian the higher order 
cumulants C q ^ 2 (xi, ■■■! x n, j + ] vanish, and we have only 

({x(x) ,x( x ')}) = C>2 (xi,X2',j + ], which is proportional 
to the Hadamard function. We shall assume that we have 
shifted the noise to a new noise with zero mean, and corre- 
spondingly we write the first cumulant as a sepa rate dissi- 
pative term. The stochastic effect ive a ction (4.2) gives the 
Langevin equations according to (3.53) as 



SS A [z] 



dx f dx' 5 4^-G R {x,x')j{z,x') 

Sj{z,x) 
oz£{t) 



(4.4) 



We have used the fact that z~ = implies j~ = 0, j + = j, 
5j~ / Sz~ = Sj/6z, and z + — z. For a system with N parti- 
cles (4.6) are N x D coupled stochastic integrodifferential 
equations for the stochastic particle worldline coordinates 
z£[t). Each realization of the stochastic field x( x ) gives a 



tions are stochastic functionals of x{ x )- I n linear models, 
the current j is proportional to the particle degrees of free- 
dom (here z), and therefore Sj/Sz^ = g u (x) is i ndependent 
of z. From this it immediately follows that (4.4) becomes a 
linear stochastic differential equatio n w ith additive noise. 
For such linear models, the mean of (44) gives the classical 
limit 



+ 2 dxdx'g nfl (x)G R (x,x')j{z,x') = 0, (4.5) 



SS A [z] 
5zit(r) 

making use of (%) = 



B. Nonlinear Coupling 

For nonlinearly coupled theories the second cumulant 
(under the Gaussian noise approximation) will depend on 
the self-consistent solution for the semiclassical equations of 
motion including the backreaction effect represented by the 
first cumulantQ. The linearized self-consistent approxima- 
tion requires expanding about the semiclassical solutions 
zfi. We define z n = z n — 
motion are therefore 



The linearized equations of 



= E/ dT mZ" m {T m ){-^ 



S 2 S A [z] 



2 / dx dx 



, 5 2 j{z,x) 
5z"Sz% 



5z^ (r m ) 8zn 
G R (x,x')j(z,x') 



+ 2 j dxil.r 
+ / dx 



,8 j(z,x) nR[ ,^ Sj(z,x) 



S -%^x(x)} + 0(z>) 



(4.6) 



The first term on the RHS is a kinetic term for the z n 
(there may also be linearized force term if Sa includes an 
external potential V (z n ). The second and third terms give 
history dependent nonlocal dissipation and multiparticle 
interactions. For brevity, we will take repeated particle 
indices (e.g., z m (5/5z m )) to always imply both summation 
over m and the integration J dT m . The fourth term is an 
additive noise. (Note that the noise is only additive for 
nonlinear theories at this lowest O (z) order.) 

These Langevin equations show the intimate connec- 
tion between particles and fields within the stochastic- 
semiclassical limit. For each realization of the stochastic 
field x{ x )i the Langevin equation has solutions z n [x, r) that 



8 For linear theories, the second cumulant is independent of the 
first, and hence linear theories do not have noise that is history 
dependent. 



16 



are functionals of x( x )- I n this sense, the quantum statis- 
tics of the field is carried by the particle trajectories. We 
could as well have treated the field as the system, and inte- 
grated out the particles, to find a stochastic effective action 
S s toch[f ' ], and Langevin equations 



SS s toch [ l ^ ± ] 



Sip 



0, 



(4.7) 



•p-=0 



where the stochastic field solutions are functionals of 
stochastic particle trajectories. In this way, the field would 
encode the quantum statistics of the particles. 

If the particle- field coupling ha d been linear, and j(x) 



had factored as /i„ M (x)z^(r), (4.6) would have reduced to 



the linear stochastic differential equation 



2 J dx J dx' h nfi {r, x)D R (x,x') 



SS A [z] 
Szft(r) 

x h mv {x')z v m [T l ) + J dx h nfl (x)x(x) 
SSa[z 



(4. 



g^,, +X nll {T) + I dr' A„ m ,^(r,r'X n (T') = 0, 



where 



$n/x( T ) = / dx h nfJ ,(x)x(x) 



(4.9) 



and 



A„ m , Ml/ (r, t') = 2 J dx J dx' h nfJi (T,x)D R (x,x r ) 

x/w(r',x'). (4.10) 



Eq. ( |4.SD are Langevin equations with additive noise 
X nfl (r) and linear, though still nonlocal, dissipation kernel 
A„ m/i „(T, t'). The equation is a general form of the many- 
particle Langevin equations for Brownian motion with non- 
local interactions and a linearly dissipative environment. 
It is equivalent to the linear N-particle detector Langevin 
equations found by Raval, Hu and Anglin || . 



C. Higher-Order Noise Cumulants 

Rather than working with stochastic field variables that 
couple to j~[z;x), it is convenient to define local stochas- 
tic forces r]% (r„) that couple directly to the coordinates 
z£ (t„) . These are just the stochastic forces that appeared 
in section III. We shall now be more explicit in our deriva- 
tion. Using the generating functional 

W[j ± (z ± )] = -ihlnS IF [j ± (z±)], 

we define the qth order cumulant (for the noise rj) to be 

a;Z+ 1 ...,Z + ] (4.11) 

s q s IF [j+,r] 



nM = r M (t 

~ °g(n...m)ft...iA''»> 



SZn {T n )...SZm{T m ) 



The probability distribution P[rj; z + ] is defined by requiring 



F\j ± } = J Dr£P[r£;z+] 



(4.12) 



x exp 
= exp 



- / drz£ (r)?? nAl (T) 



n 



z^~ z v ~ 



X < ^2('nm)/ii/(' r ™' Tm ' Z + ] + •■• 

+ E [dr n ---dr p z 



x C 



(v) 



n...p 



(t„,...,t p ;z + ] + ....}. 



P[rj; z + ] is the characteristic functional of F [j^] with re- 
spect to zti~ ', and C, li "'\axe the cumulants of the noises 

" n 1 q{n...m) 

rj n ^{T n ). Note that there is a noise for each index // (in- 
cluding the timelike component), and for each particle n. 
Because the influence functional is not Gaussian in z, the 
higher order cumulants C q y 2 do not vanish. The mean 

(?]nfi) of the noise is the first cumulant Cjj^ {r n ', z + ]. It 
gives the dissipative term, and as before, we will shift the 
noise so that rj has zero mean. Then the stochastic effective 
action becomes 



5^,77] = S^} I drz^(T n ) (4.13) 

n J 

"WS),(M + ]+^(r)). 

The stochastic equations of motion follow immediately 
from (6S v /5z£-) \ z - =0 = 0- They are 



(4.14) 



If (4.14) were linear, we would interpret it as a Langevin 
equation for the particle trajectories with additive, but col- 
ored noise i](t). In fact, as it stands, (4.14) is nonlinear, 
and the noise depends in a complicated way on the trajec- 
tories z because the probability distribution P[r), z + ] is a 
functional of the trajectories. As before, we will assume 
that decoherence suppresses large fluctua tions around the 
mean-trajectories, allowing us to expand (4.14) in terms of 
the fluctuations z around the solutions z defined by 



SS z [z] 



(4.15) 



;=o 



We find the linear stochastic integrodifferential equation 
for the fluctuations z(t) — z(t) — z(t) : 
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dr„ 



6 2 S z [z] 



E 



^ / "' m U:;;. T M:;;,(^! 

(»?) 



5Cy> (t,z\ 



where we have used ( [4.15 ) to eliminate the non-fluctuating 
terms. If we let the indices denote the set of indices 

(n, (i, t') , with all integrations a nd su mmations implied by 
repeated indices, we may write (4.16) as 



where 



and 



(Tij + Rij)zj + r)i 



T i: j = (S 2 S z [z]/5ziSzj 



Rij - 



(sc[f[z]/Sz 



(4.17) 



(4.18) 



(4.19) 



is a kinetic term, plus a possible external potential 
term. The diagonal terms of Rij represent dissipative 
forces, while the off-diagonal terms are nonlocal particle- 
particle interactions, called propagation terms in M fo r the 
case of the linear particle-detector model. In (4.16), we 



have made a Gaussian approximation to the noise by ne- 
glecting all higher order cumulants and evaluating the sec- 
ond order cumulant Crp (n , T2 ; z + ] on the mean-trajectory 
z. The noise correlator (matrix) for T)i is then 



(4.20) 



The diagonal terms of the noise matrix gives single particle 
noise, of the general type found in QBM models. The off- 
diagonal terms give what Raval and Hu have called corre- 
lation noise, representing the nonlocal correlations between 
different particles as mediated by the field, despite it being 
coarse-grained away. 

Finally, we should show that the description in terms 
of stochastic fields x is equivalent to that in terms of lo- 
cal stochastic forces rj. Making use of the form of SifU^], 
we may immediately evaluate the first cumulant, which is 
given by 



Ci„(r, z] = dx dx 



j Sj{x,z] 
Sz^(t) 



G R (x,x')j(x',z]. (4.21) 



This gives the same result for dissipation as was found in 
Eq. (4.4). The noise r\ is 



dx 



Sj(z,x) 
5z£(t) 



(4.22) 



We find the correlation of 77 by direct computation: 



(^(rMr')) 



dx dx 



8z^-(t)8z v -(t') 

,f Sj(x,z) Sj(x,z] 
L Sz^(t) 5z u {t') 
Sj(x,z]6j(x,z] ^ H 
6z*{?) 5z»{t) J 1 ' ' 
,5j{x,z] Sj{x,z) 



(4.23) 



dx dx 



5z^( T ) Sz u (t') 



(x(x)x(x')) , 



where we ha ve used the symme try o f the Hadamard func- 
tion in ( 4.23 ). This agrees with ( 4.22 ). Therefore, one finds 
the same stochastic equations of motion whether one works 
in terms of a stochastic field x> extended though space, or 
a stochastic force 77 acting solely on the particle worldline. 



V. GENERALIZED FLUCTUATION-DISSIPATION 
RELATIONS 

Fluctuation-dissipation relations (FDRs) play an impor- 
tant role in both equilibrium and nonequilibrium statistical 
physics. They are vital in demonstrating the balance nec- 
essary for statistical dynamical stability of a system in the 
presences of external fluctuation forces. Because of the in- 
timate connection between statistical, quantum, and rela- 
tivistic processes in which fluctuations and dissipation play 
a role, we now address the question of whether fluctuation- 
dissipation relations exist for the more general nonlinear 
particle-field models under study. This inquiry will be im- 
portant for understanding the stochastic limit of particles 
in quantum fields, and the relationship between dissipation, 
stochasticity and vacuum fluctuations. 

For linear quantum Brownian motion, the Langevin 
equation takes the form 



dt'fj,(t,t')q(t r 



+ cu 2 q 



(5.1) 



The fluctuation-dissipation relation is an integrodifferential 
identity relating the dissipation kernel /j,(t, t') and the noise 
kernel v(t,t') = (£(t)f (*')),, b y a universal kernel K(t,t') 
whose properties are independent of the spectral density 
of the environment and the coupling strength e. For QBM, 
defining the kernel j(s) by /i(s) = d r y(s)/ds,the fluctuation- 
dissipation relation is 



ds'K(s,s')j(s) 



(5.2) 



where 

K(s,s') 



kdk ( 0kh\ 

coth — — cosfc (s — s ) . (5.3) 

7T V 2 I 



When the dissipation is approximately local, and the noise 
approximately white, we recover the familiar Einstein- 
Kubo relation W2[ . A special feature of the ordinary FDR 
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for this condition is that the kernel K is independent of 
the system variables. This turns out not to be true for 
nonlinear theories. 



For the full nonlinear Langevin equation (4.14) found 
in the previous section, any FDR's that are found should 
take the form of a relationship between the cumulant 
n(L( T " ,2 l' gi ym g radiation reaction, and the noise cor- 
relator for rj ntl (t„) . But for the nonlinear Langevin equa- 
tion this connection between noise and radiation reaction 
is much more complex than the familiar one characterizing 
the linear response regime. In paper II, we see an example 
of the cou nter- intuitive behavior of the nonlinear Langevin 
equation (4.14). There it is shown that for uniformly ac- 
celerated particles, the radiation reaction term given by 
the first cumulant (which is often conveniently yet wrongly 
thought of as dissipation) vanishes despite the presence of 
quantum fluctuation-induced noise. The relationship be- 
tween radiation reaction and the noise embodied in any 
kind of general FDR is of a subtle nature, and it is not 
correct to conclude, on the basis of one's experience with 
the linear response or equilibrium regimes, that radiation 
reaction is flatly balanced by fluctuations. What is bal- 
anced by fluctuations are the variations in the radiation 
reaction force about the average (i.e. classical) force asso- 
ciated with deviations of the particle trajectories around 
the average (i.e. semiclassical) trajectory. We identify the 
stochastically varying forces associated with fluctuations in 
the trajectories as the true dissipative effect, and call this 
quantum dissipation because of its quantum origin (i.e. this 
effect vanishes when h — > 0). Therefore, the proper state- 
ment is that quantum dissipation is balanced by quantum 
fluctuation induced noise. 

To see the emergence of the usual kind of FDR within the 
linear response regi me, w e should consider the linearized 
Langevin equation ( 4.17 ) for fluctuations z M around the 
semiclassical solution z^. Such a FD R wi ll be a rela tion 
between the matrices Rij and Cy in (4.19) and (4.20), re- 
spectively. However, we can immediately see that there 
will not be a completely general FDR involving all compo- 
nents of these matrices B. The matrix Rij is dependent 
on the field commutator evaluated at the spacetime points 
Xi and Xj. For spacelike separated points, the commutator 
identically vanishes. On the other hand, the matrix C'ij 
is proportional to the anticommutator evaluated at these 
points, which does not vanish. For example, two oppositely 
charged particles in a uniform electric field sufficiently sepa- 
rated will never enter each other's causal future, and there- 
fore R12 will always vanish. However, the particles will be 
correlated through the nonlocal vacuum fluctuations of the 
quantum field. These nonlocal correlations between the 
two particles will be determined by a non-vanishing N12 
(see Figure 2). 

In exploring the question of what kind of FDR relations 
do obtain, we generalize the treatment in RHA j^] for par- 
ticle detectors on fixed trajectories linearly coupled to a 
quantum field. We will find that an FDR looking superfi- 
cially like the one found in M does obtain, but that it is in 



fact a much more complicated relation. We first define the 
retarded matrix Gf- = G R (zt,Zj) . This matrix is the same 
as flij = Re(Zy) in || except for the important difference 
that here the matrix is a nonlinear function of the particle 
variables, whereas in RHA, the kernel is independent 
of the particles dynamical variables. 

We first consider the diagonal elements of G R , which, 
in part, are responsible for radiation-reaction forces on the 
ith particle. When the field is massless, radiation from the 
particle will travel outward on its future lightcone. Assum- 
ing that the particle trajectory Zi is everywhere timelike, it 
will then only interact with its own field at the instant of 
emission. It is straightforward to see that as a consequence 
of this, Gf[ oc S(ri — t[), and hence is local. Actually, these 
conclusion depend on the spacetime dimension. In 2+1 di- 
mensions, the massless, retarded Green's function is not 
restricted to the lightcone, and so radiation reaction is not 
local. In 1+1 dimensions, infrared divergences are a prob- 
lem because the retarded Green's function does not fall-off 
with distance, and therefore, radiation fields (and radiation 
reaction) must be carefully defined (see p3[). 

In 3+1 dimensions (unlike 1+1), we must confront an 
additional complication: G R diverges as t- — > t;, and thus 
does not give a finite radiation-reaction force as it stands. 
We follow the approach of writing the kernel Gf- as 



1 



1 



and define the new matrices 

f 

2 

and 



= 7;l (j y — ^ij) 



1 



(5.4) 



(5.5) 



(5.6) 



Gfj is the advanced Green's function evaluated at i and 
j. The matrix Gfj is given by the field commutator 
G c (x',x) = ([ip(x'), <f(x)]} evaluated at x' = Zi and 
x = Zj. If we define the radiation fields by ip rad (tf) = 
<p out (t f ) - <p m (U), then 



<p rad (x) = / G c (x,x')j(x')dx'. 



(5.7) 



The energy carried to infinity by ip rad is consistent with the 
Larmor formula (or its scalar field equivalent). But while 
the Larmor formula gives the net radiant energy assum- 
ing that the particle is un-accelerated at initial and final 
times (or that the particle is undergoing periodic motion), 
it does not account for all instantaneous forces acting on 
the particles. In particular, there are both short-range po- 
larization type forces (including the so-called acceleration 
or Shott fields) acting between particles and modifying the 
radiation reaction forces giving corrections to what would 
be expected on the basis of the Larmor formula alone B. 
Neither of these additional types of forces do net work, but 
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they do modify the equations of motion. Since we are inter- 
ested in the detailed, finite-time evolution of the particles, 
these type of forces should not be neglected. 

In fact, only the off-diagonal elements of the matrix £r£y 
give Coulomb- like (or polarization) effects. The local, but 
divergent diagonal elements Gf i: renormalize the particle 
masses. Both of these effects are independent of the quan- 
tum state of the field, since they depend on the field com- 
mutator, and hence do not require a quantum mechanical 
treatment to understand. 

In Paper II , we will discuss renormalization effects in 
detail; here, we assume that renormalization has been per- 
formed in order to focus on the FDR question. We assume 
that the diagonal G? dependent terms have been combined 
with the kinetic terms to give mass renormalization, and 
we define the dissipation matrix 



A.. — r< H _ c F X- . 



(5.8) 



Similarly, we define the noise matrix 



Ni 



G H (zi, Zj) 



which is found from the field anticommutator. 
Restricting to the single (i th ) particle, we define 



dr 



zfd],A l 



(5.9) 



in analogy with the usual FDR in QBM. The chain rule 
has been used to write d/dr = zf(dldzf) = zfd],. In 
the expressions ( |5.8|) and Q5.9| ), the index i is not summed 
over, and the derivative is only applied to the first z, term 
in An (t, t') = A[zi (r) , Zi (r')). Our definition and place- 
ment of the derivative agrees with , where one can find a 
discussion on why the appropriate placement of the deriva- 
tive depends on whether the interaction term is written as 
derivative coupling, or minimal coupling. Our case is a 
modified example of derivative coupling. 

We may now demonstrate the existence of a FDR for 
the diagonal elements of An and Nn by making use of the 
general properties of the commutator and anticommutator 
for a field in an initial thermal state at temperature T. For 
the diagonal elements, we have An (t, t') = Gf { . We will 
now drop the i index altogether, writing An = A, Nu = N, 
and Zi = Z, since we are concerned with only one particle 
in the following. Defining x^ (r, r') = z^(r) — z^(t') and 
K = |k| , and using the Fourier transforms of the various 
Green's functions, we have 



(5.10) 



N{t,t') = J dke- tk » x 5(k 2 ) 

x (9(k°) - 9(-k°) + 2n(«)) 



dke 1 



2n(n) 



where n (k) is the Bosc distribution term 

= r Tu r\ — T- 

exp [k/kb-l \ — 1 

A (t, t') is similarly given by 



A(t,t') = / dke lk x 



-\-%KX 



2 k 



(5.11) 



(5.12) 



without the thermal Bose distribution factor. We assume 
that the particle trajectories are timelike so that x is time- 
like. Then, after taking the gradient of A, we may choose a 
coordinate system so that x = 0. We find the simple result 



V x ^ 3+1 (z)|x=o = Jdkk 



( sinrexo 



V K 

since the integral is odd in k, while 







(5.13) 



_d_ 

dx 



-A 3+1 (x)\ x=0 = 47T / dKK 2 COs(KX ) 

o J 



= -An 



dx 2 



dn cos (kxq) 



o ./o 

d 2 



= -l^ 2 ^S(x ) (5.14) 
ox 

in 3 + 1 dimensions. For comparison with other work, we 
also note that in 1 + 1 dimensions d/dxo (A 1+1 (x)| x= o) = 
4:Tt5(x ) and V X A 1+1 | X=0 = 0. Use of 



S(z°(r) - z°(r')) - S(r - r')/ \dz°/dr\ T=T , 



(5.15) 



requires the trajectories to be timelike, which we assume is 
the case, at least for z (r) . Therefore, z° (r) is a monotonic 
function of t, z (t) > 0, and we may use 



1 



d 



1 d 



dz% zo (r) dr \ z (r) dr 



(5.16) 



In then follows that (5.9) is given by 



T 3+1 = 16^e 2 i ^(z°(T) - z°(t')) 
dz 2 



= \Qn 2 e 2 ^- 

OT 



1 d (I 

Zq OT \ Z u 



in 3 + 1 , and by 

r 1+1 = -4ne 2 z 6(z (T) - z°{r')) 
= —4ire 2 5(T — r'). 



(5.17) 



(5.18) 
(5.19) 



2 k 



in 1 + 1. This result holds for any state of the field since it is 
derived from the field commutator. Since d^A is covariant, 
this result is in fact coordinate independent, despite our 
use of x = in the intermediate steps of the derivation. 
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For the noise matrix element N, we have 



N = I dn 



;6(l + 2n(«))|e 



-i«(x°-|x|) + e+4K (x + |x|)| ^ 

(5.20) 

where b = Airn in 3 + 1, and b = in 1 + 1 dimensions, 
respectively. Finally, we define the fluctuation-dissipation 
kernel 

f°° dn 

K(t,s) = / — &[l + 2n(«)] 

x -[cosk (v(t) — v(s)) 

+ cosk(u(t) - u(s))} (5-21) 

in terms of the out-going coordinates u = x — |x| and the 
in-coming coordinates v = x° + |x|. We may now verify the 
fluctuation dissipation relation, 

/>oo 

N[z;t 7 t')= dsK[z;T,s)T[z;s,T'), (5.22) 



by substituting K, and using the locality of T. Eq. ( 5.22 ) 
f ollow s immediately by virtue of the delta function in Eq. 
(5.18) for r 1+1 . In 3 + 1 dimensions, we integrate by parts 
thereb y sh ifting the operator d 2 /dz 2 from r 3+1 to K in 
Eq. (j5.22| ). The proof of the FDR then follows imme- 
diately, once again, by virtue of the delta function. Kn 
agrees with the kernel found by Raval, Hu and Anglin for 
T = in 1 + 1 dimensions This FDR generalizes their 
result to the more physical 3 + 1 dimensions, to non-zero 
temperatures for the quantum field tp, and to the case of 
nonlinear interactions. 

Of particular importance is that unlike the traditional 
FDR (which is a universal kernel independent of the system 
dynamics), the kernel K here is a functional of the system 
history Zj. Hence, (5.22) generally represents a complicated 
(nonlinear) relationship between the noise and dissipation 
matrices. We have alrea dy commented on the significant 
difference between ( 5.22| ) and the ordinary FDR that ob- 
tains for the linear response regime, one example being the 
vanishing of the (average) radiation reaction force for uni- 
formly accelerated charges despite the existence of quan- 
tum fluctuations. This indi cates that the formal relation- 
ship as embodied by ( |5.22 ), while looking like that kind 
of FDR found for linear theories, must be carefully inter- 
preted. 



corresponding quantum operators for the particle space- 
time coordinates, when the coarse-grained particle trajec- 
tories are sufficiently decohered. Specifically, the equations 
of motion follow from the generating functional for world- 
line correlation functions. When decoherence is sufficiently 
strong, quantum features associated with higher-orders 
quantum fluctuations (the (D(z~ ) terms in the CGEA) are 
exponentially suppressed. After doing the z path integrals 
to one-loop, the remaining z + path-integrations are equiv- 
alent to a stochastic path integral with respect to some 
system-history-dependent non-Gaussian probability distri- 
bution. Hence, one finds that the quantum dynamics of 
correlation functions are approximated by non-Markovian 
stochastic dynamics with colored noise- the quality of the 
approximation depending on the degree of decoherence. 

The resulting nonlinear Langevin equations for the par- 
ticle spacetime trajectories may be expressed in a num- 
ber of equivalent ways: in terms of a stochastic force rj (r) 
coupling directly to the particle coordinate, a stochastic 
variable 2 coupling to the functional f[z,t), or a stochas- 
tic field x ( x ) coupling to the conserved current j (x) . A 
cumulant expansion may be used to systematically incor- 
porate higher-order noise contributions. At order e 2 (low- 
est order in the cumulant expansion after the dissipation 
term) the noise is Gaussian. At high temperatures, the 
approximately Gaussian noise is effectively white. As one 
moves from the classical to stochastic regimes, particularly 
at low temperature and for short-times, higher order (non- 
Gaussian) noise plays an increasingly important role. Mov- 
ing deeper into the quantum regime (when decoherence is 
less complete), higher-order quantum corrections become 
significant and the stochastic approximation outlined here 
fails. 

In Paper II, we will apply the formalism developed here 
to the specific problem of radiation reaction for spinless 
particles. We derive a Langevin equation whose average 
trajectory is well-approximated by the Abraham-Lorentz- 
Dirac equation at late times, but where time-dependent 
renormalization (arising from the nonequilibrium quan- 
tum treatment) modifies the early time behavior making 
it pathology free. 
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VI. CONCLUSION 

In this paper we have presented a general framework 
for describing stochastic particle motion in quantum fields. 
These classical nonlinear stochastic equations of motion are 
good approximations to the underlying quantum theory, in 
the sense that the noise average of the stochastic variables 
are nearly identical in form to the expectation values of the 



APPENDIX A: INITIAL- VALUE METHODS 

The Closed-Time-Path (CTP), Schwinger-Keldysh, or 
'in-in' method P,pj,|6|,|l4j is designed for the computation of 
the evolution equations for correlation functions of quan- 
tum operators. This is in contradistinction to the compu- 
tation of 'in-out' transition amplitudes as in scattering pro- 
cesses, or in dealing with equilibrium conditions, where one 
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either knows about, or knows what to expect in, the final 
state. As such, it is particularly appropriate for nonequi- 
librium processes in quantum field theory, characterized 
by [pin,H] 7^ 0, and for which the final state is generally 
unknown. Even if one knows or assumes that the system 
does eventually come to equilibrium, how it approaches 
equilibrium requires some knowledge of the evolution dy- 
namics involved. Because the CTP formalism explicitly 
works with true expectation values (in\0 \in) of a Heisen- 
berg operator evaluated with respect to initial states, it 
produces equations of motion that are both real and causal 
EJ. These are essential properties that are lacking in the 
equations of motion found from the "in-out" effective ac- 
tion. When there is a natural division between system and 
environment (relevant and irrelevant) degrees of freedom, 
the influence functional invented by Feynman and Vernon 
]Tl[ ] is most appropriate as it explicitly implements a coarse- 
graining of the environment. It has been shown to be inti- 
mately related to the CTP coarse-grained effective action 
method (see Appendix C). 

The influence functional is useful when we treat a dis- 
tinguishable system (corresponding to the relevant degrees 
of freedom) in contact with an environment (the irrelevant 
degrees of freedom). Such situations arise frequently in 
practice, though from the point of view of the microscopic 
degrees of freedom, the division is not always so clear-cut. 
Indeed, we never directly observe the most elementary de- 
grees of freedom, and thus it is important to find the effec- 
tive or collective variables which capture the most relevant 
physics of the system at the observation scale in question. 
In this work, the mass difference between the particles (e.g. 
electrons) and the massless field quanta (e.g. photons) pro- 
vides the basis for a reasonable partition between system 
and environment. For the problems under study here, we 
will choose to measure only particle observables, leaving 
the field degrees of freedom as the environment; the physi- 
cal appropriateness of this choice of the partition depends 
on the particular parameters of the theory and the physical 
questions one is asking. In other contexts it may be more 
natural to measure field observables, treating the particle 
degrees of freedom as an environment [ fl5| . 

Our use of the influence functional in this work involves 
the evolution of initially uncorrelated states. Despite the 
mathematical convenience of this assumption, we must 
identify the physical conditions under which this ideal- 
ization is appropriate. In fact, a sharp system and envi- 
ronment division of a particle interacting with a field is 
unattainable: such factorized states do not exist in infi- 
nite dimensional Hilbert spaces since it may be shown that 
infinite energy is required to produce them from the true 
(interacting) ground states ||. A field cutoff makes fac- 
torized states well-defined mathematically, and reveals the 
need for regulating the nonlocal particle-field correlations 
in the theory in addition to the local renormalization effects 
that arises in more traditional field theory applications. 

Time-dependent renormalization is a generic feature of 
finite-time nonequilibrium processes, because the system 



may not have enough time to equilibrate with the environ- 
ment. Operating at a faster scale in the dynamics will be 
the time-dependent dressing of the particle j4(|. It is not 
sufficient to renormalize the parameters of the theory in 
the usual way, because the renormalization procedure only 
provides local counterterms. Here, we are explicitly inter- 
ested in the nonlocal correlations. For example, even with 
the natural division of system and environment variables 
given by the mass gap between the massive particle and 
a massless field, the dressed system still involves extensive 
nonlocal correlations in addition to the more familiar local 
shift in the particle mass. Furthermore, at high energies, 
the energy gap between massive particles and massless field 
quanta is unimportant, and the distinction between parti- 
cle and field variables becomes blurred. 



1. Path integral derivation of evolution kernels 

The system variables z£ are the spacetime coordinates 
of the worldlines of a collection of charged particles (the 
subscript denotes the n th particle), each independently 
parametrized by r„ . The environment /field tp (x) is equiv- 
alent to a collection of harmonic oscillators. We employ 
a field description for the environment of massless parti- 
cles because radiation is a significant feature, even at low 
energy, due to the zero-mass threshold for the production 
of massless field quanta. On the other hand, at low ener- 
gies, creation of massive charged particles is strongly sup- 
pressed. Therefore, we opt for the worldline description of 
massive-charged particles to take advantage of its particle- 
like features. The description of particle creation, and the 
handling of states with indefinite particle number (such as 
coherent field states), are conveniently depicted by field- 
theory, but particle creation, and related processes, can 
also be addressed in the worldline framework j4?J]. The 
naturalness of the worldline framework can be seen from 
a background field perspective, where the mean particle 
trajectory is analogous to a mean-field. 

We call the closed system of particles plus a field the 
universe, and its action will be assumed to be of the form 

S[z, ip] = S z [z] + S v [p] + S int [z, ip}. (Al) 

In most studies of quantum open systems using the influ- 
ence functional the interaction term Si n t is linear in z and 
ip, while the terms S z and S v , determining the isolated dy- 
namics of the two sub-systems, are quadratic [ppp[P|. A se- 
vere limitation of quadratic theories with linear coupling is 
that they are equivalent to non-interacting free theories in 
terms of appropriately redefined system and environment 



9 Hu, Paz, and Zhang have considered in detail quantum Brow- 
nian motion (QBM) with nonlinear (polynomial) particle-bath 
interactions using the influence functional and perturbation the- 
ory @. 
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variables that diagonalize the Hamiltonian. Such models 
do not explain how a decoherent basis can arise dynami- 
cally, and hence are limited in their ability to account for 
emergent classical and stochastic behavior from the effect 
of noise-induced decoherence. 

We work in a initial value formulation that is appropri- 
ate to the study of nonequilibrium quantum open systems, 
with the aim of deriving the propagator that gives the evo- 
lution of the reduced density matrix p r ^ for a system ini- 
tially prepared at time ij. (In the following, the subscripts 
i/f will denote initial/final.) More generally, we should 
consider the evolution of p r .T.t initially defined on some ini- 
tial Cauchy hypersurface Sj. The examples in this paper 
involve flat (Minkowski) space, but the methods are imme- 
diately generalizable to curved spaces |m|,[37| ■ The reduced 
density matrix associated with coarse-graining the environ- 
ment field allows us to evaluate the system observables at 
spacetime points after E^. To derive its propagator we be- 
gin by considering the full evolution operator for the closed 
system. We assume (as part of the initial condition) that 
the particles begin on the surface Sj with parameter val- 
ues Tj, where z° (r,) = U. The initial state of particles plus 
field is therefore defined on Sj. 

We assume a unitary operator Ufi giving the evolution 
of the initial density matrix pu,iOi the closed universe ac- 
cording to 



where the measures are defined as 



puj = Ufipu,iU f fi = JfiPu 



(A2) 



where Jfi, the density matrix evolution operator, is defined 
by the above expression. We take \z) to be the (spacetime 
position) eigenstates of a Schrodingcr picture operator 



such that 



(we drop the particle index n 



for brevity). The operator z M is conjugate to the momen- 
tum operator p^, whose eigenstates are momentum states 
\p) . These operators satisfy the "equal time" commutation 
relation [z^p 1 '] = ihg^ v . Momentum states are physical 
if they are on-shell: p^p^ = m 2 . The spacetime position 
states \z) are not in the physical Hubert space because they 
do not satisfy the Hamiltonian constraint associated with 
the reparametrization-invariant (i.e. generally covariant) 
relativistic particle. However, they do have overlap with 
the physical Hilbert space, and therefore may be used as 
a convenient basis despite not depicting physical particle 
states. 

Similarly, let \<p) be the eigenstates of the Schrodinger 
operators ip = <p(ti) such that tp \tp) = ip \tp) . We may pick 
as a basis for the universe's Hilbert space Hu = 7~Cz <S> T~t v 
the direct product states \zp) — \z) <£> \p) . We now define 
the transition amplitude (i.e. the matrix elements of Ufi) 
to go from the initial state \zipi) to the final state \zfipj) 
by the path integral expression 

r z ffif ( 1 

(zf<pf\Ufi\zi(pi) = J DzDpexp ^-S[z,p] 



K(zf,(pf,Zi,ipi), 



(A3) 



J x,t J -°° 

Dz = n / dz ° M / dz ( 1 



(A4) 
(A5) 



Actually, this discussion of the quantum theory leaves 
out many important details because the relativistic par- 
ticle has a reparametrization-invariant action, and is 
therefore a gauge theory. The particulars of how the 
reparametrization-invariant theory is quantized depends on 
the form of the action S z , an d how the gauge redundancy 



is handled. As it stands (A3) is ill-defined since it involves 
summing over many gauge-equivalent worldlincs. Since we 
are interested here in the semiclassical particle behavior, 
we will make the simplifying assumption that the world- 
line path integrals have been suitably gauge-fixed. If there 
is residual gauge-freedom (such as the lapse N that ap- 
pears in the proper-time path integral representation) than 
the measure for these variables will be implicitly included 
in the definition of Dz. In applying the stationary phase 
method to define the semiclassical solution, we will then 
fix any residual gauge-variables to the "classical" (i.e. ex- 
tremal) value. For instance, the proper-time path integral 
involves an integ rati on not only over worldlines as given by 
the measure in ( |A5| ) , but also an integration over all pos- 
sible proper-times for a particle to go from Zi to Zf. How- 
ever, the stationary phase approximation is dominated by 
the worldline whose parameter is the classical proper-time 
satisfying 



(dz^ (r) /dr) 



1. 



(A6) 



For simplicity, we shall fix the parametrization to the clas- 
sical proper-time value so that the semiclassical worldline 
solution satisfies ([A6|), and wc shall ignore the integra- 
tion over residual gauge freedom. Doing this therefore ne- 
glects quantum fluctuations in the value of the lapse N. 
Since such fluctuations may also have a stochastic charac- 
ter that manifests in the stochastic regime, we should in- 
clude them in our consideration for a complete description 
of the stochastic regime. 

There are additional subtleties to this worldline path 
integral regarding the class of worldlines that are consid- 
ered (this is related to the form of the action S z and the 
choice of gauge). In our second series, concerned with 
fully quantum relativistic effects, we employ the so-called 
proper-time gauge. In that gauge, one includes all possi- 
ble worldlines going between the spacetime points Zj and 
Zf, including spacelike and backward-in-(Minkowski)-time 
worldlines. We will discuss then why, and how, a restricted 
path integral must be employed where worldline paths go- 
ing back in time before t% must be excluded, but all other 
possible paths are included. This restriction is given in the 
path integral measure J Dz' 1 where the integration range 
of z° is ti < z . We mentioned earlier how this temporal 
restriction leads to a kind of finite-size effect modifying the 
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particle dynamics near the initial time hypersurface. It is a 
nice feature of the spacetime (worldline) formulation that 
finite-size (e.g. boundary) effects coming from both spa- 
tial and temporal restrictions of the path integral may be 
incorporated on equal footing. 

We may alternatively employ the worldline analog of the 
Coulomb gauge where t = r. In this gauge, particle trajec- 
tories go strictly forward in time t, and we can set t = r n , 
for all n. Such a gauge choice, which is incquivalent to 
the proper-time gauge choice, produces the Newton- Wigner 
propagator for relativistic particles, rather than the Feyn- 
man propagator. When particle creation processes are im- 
portant, the difference is significant. Here, we assume a 
priori that we are in a regime where particle-creation effects 
are insignificant. In the Coulomb gauge (t = t), z° is no 
longer an independent degree of freedom, and the Langevin 
equations derived in the body of this paper should be mod- 
ified by z M — > z % (and r — > t). 

In fact, insofar as this paper is concerned, we could as 
well have adopted the semiclassical paradigm of treating 
the particle worldlines as classical variables. We could then 
just take the coarse-grained effective action found by inte- 
grating out the quantum field degrees of freedom as the 
definition of a semiclassical (or stochastic) theory. If we 
had followed this philosophy, the results of this paper would 
have been the same, but we would have lost precious in- 
sight into the meaning (and range of applicability) of the 
semiclassical or stochastic theory. 

Writing the initial density matrix as in (|2.2| ), the final 
density matrix is then defined by 

Puj(zf,z' f ,tp f ,tp' f ) (A7) 

x J{zf,z' f ,ipf,<p'f;zi,z , i ,<pi,ip' i ) 
x p u {z t ,z' l ,ip, l ,ip' l ). 

The density matrix evolution operator is given by 

J(zf,z' f ,<Pf, i P'f^u^'Ph l P'i) (A8) 
= K{zf, <pf,Zi, ipi) K*(z' f ,<p' f , z'i, if',) 

= I I DzDtpDz'Dtp' 



xexp^ ~(S[z,<p] - S[z',<p']) 



2. Reduced density matrix 

Working with the density matrix evolution operator al- 
lows one to treat the full quantum theory in an initial value 
framework. However, in many contexts the full quantum 
theory is not what one is really interested in. Realistic 
situations involve the measurement of a subset of the uni- 
verse's degrees of freedom. One direct way to implement 



this physical situation is to work with the reduced density 
matrix for the relevant variables (the system) by tracing 
out the unobserved variables (the environment). Work- 
ing with the reduced density matrix implements one type 
of coarse-graining; a discussion of more general types of 
coarse-graining may be found in ]Tg| . 

Our reduced density matrix formalism follows the effec- 
tive dynamics of coarse-grained histories whose members 
are fine-grained particle trajectories but completely coarse- 
grained field histories^. Through the influence functional 
one may construct an effective action Scgea for these 
coarse-grained histories [jl2). It is important to keep in 
mind that the coarse-graining of environment variables 
alone is not sufficient to allow the description of system 
histories by classical dynamics - additional coarse-graining 
(smearing of the trajectories) is required to find truly de- 
coherent sets of histories . It is through coarse-graining 
that the statistical aspects characterized by decoherence, 
dissipation, and noise arise. 

The path integral approach gives a convenient represen- 
tation for the evolution kernel of the reduced density ma- 
trix J r found by tracing over the final time environment 
variables in (A7). One obtains 



p r j{z f ,z' f ) 

dl PfPuj(zf,z' f ,ipf,ip f ) 



(A9) 



dzidz'i diff dipidip'i 

X J(Zf,Zf,(pf,(pf,Zi,4,'Pi,'Pi) 

= J dz t dZi J r (zf, z' f ; z i7 z-)p r (z 4 , z[). 



Equation (A9) may be taken as the definition of 
J r (zf, z'p Zi, 4), which in general depends on the initial 
state p r {zi, z[). Closed form expressions may be found for 
J r only for special initial states pu,i- The most easily 
treated examples are systems and environments that are 
initially uncorrelated, and which therefore have initially 
factorizcd density matrices 



10 The finest-grained histories are the individual histories of the 
full set of quantum variables that appear in the path integral 
for the propagator of the universe [If9| (The classification and 
enumeration of histories may also be formulated in terms of 
projection operators.) Fine-grained histories are never classical 
since they will always interfere with other fine-grained histories 
that are sufficiently nearby. A coarse-graining is a specification 
of sets of fine-grained histories that is exclusive (no fine-grained 
history belongs to more than one coarse-grained history) and ex- 
haustive (all fine-grained histories belong to some coarse-grained 
history) . 
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PU,i = Pz,i ® (><p,i ( Al 0) 

[see (^)] and evolution kernel J r given by [see ( |2.5| )] . The 
influence functional is given by F[z,z'] and the influence 
action is defined by Sif[z, z] — — ihlnF[z, z'\. The coarse- 
grained effective action is defined by 



More generally, it may be shown using the representation 
( A15| ) that \F[z, z'\ < 1 for all z. This property plays a 
crucial role in decohcrcnce, and the emergence of a semi- 
classical stochastic limit. 



APPENDIX B: QUANTUM BROWNIAN MOTION 
AND OTHER MODELS 



Scgea [z, z \ = S z [z] - S z [z'\ + S IF [z, z'\ . 



(All) 



An alternative form for F[z, z'\ may be found by consid- 
ering the evolution operator for the environment degrees 
of freedom interacting with the system variables treated as 
if they were a classical (c-number) source. Defining the 
time-dependent Hamiltonian 



HJz\=H v 



H in t [fi,z], 



(A12) 



we use the formal solution for the evolution operator: 

(H v [0\ + H int [0,z])\ . (Al:S) 



U[z] = Texp<{-- / dt [H vV 



The field evolution therefore depends on the system his- 
tory z. Writing \<p/,z) = U[tf,U;z]\(pi), and using the 
expression 



(tp f \ U[z] fa 



ft 



Dip exp {t( s A 



iSint [<P,Z]) 

(A14) 



allows us to write the influence functional in a representa- 
tion independent form as 



F[z,z']= J diff dipid(p' t (<Pf\U[z]\<pi) p v (<Pu<Pi,U) 
* (<Pi\UW\<Pf) 

= Tr v {U[z]p v (t i )tf[z'}). (A15) 

When the environment is initially in a pure state, p v (ti) = 
one obtains the simple result 

F[z,z'\ = (^\W[z'}U[z] fa) = (i>' f ,z'\i> f ,z). (A16) 

showing how the influence functional measures the inter- 
ference between the different environment final states that 
would result through interactions with different possible 
system histories z. 

There are some general properties of the influence func- 
tional (and the corresponding CGEA) that are w orth not- 
ing. They follow from the representation (A15) and the 
cyclic property of the trace: 

Scgea[z, z'] = -Scgea[z',z]* and S C gea[z,z] = 0. 

(A17) 

When the interaction between system a nd env ironment are 
turned off, it immediately follows from ( A15) that 



Si F [z,z'] = 0^> F[z,z'\ = 1. 



(A18) 



1. Time dependent linear QBM 

An important application of the influence functional 
method is linear quantum Brownian motion, for which ex- 
act results are known. We review here the case for time- 
dependent frequencies and coefficients following the treat- 
ment of Hu and Matacz 0] . These results are directly ap- 
plicable to a number of interesting examples where particle 
creation via parametric amplification is important, such 
as in quantum/atom optics and quantum field theory in 
curved spacetime, such as Hawking and Unruh radiation. 

The most general action for which the exact influence 
functional may be found is 



S[z,q] = S[z} + J2 J ds 



+ f n [z; s)q n {s)) , 



(Bl) 



with stipulated time-dependent masses m n (s), frequencies 
tu n (s) , and cross-terms b n (s) . The function f n [z;s) is a 
functional of the system variables, and may be left arbi- 
trary. When the system and environment arc initially fac- 
torizable {p zq = p z ® p q ) the influence functional is given 
by 



F[z, z'] = exp 



ds / ds' 

Jti 



(B2) 



X (fn[z;s) - fn[z']s)) 

x (( n (s, s')f n [z, s') + Q(s, s')f n [z', S')) } 

where the nonlocal influence kernel £ n (s,s') encodes the 
complete influence of the environment on the system. In 
(B2), the real (v) and imaginary (p) parts of the kernel £ 
are called the noise and dissipation kernel, respectively. 

It is significant that only the imaginary part (the noise 
kernel) is modified by the environment being at a non-zero 
temperature. The dissipation kernel modifies the extremal 
solution of the path integral for the evolution operator, and 
imparts a dissipative component in the classical equations 
of motion. The noise determines the diffusion term in the 
master equation, which is responsible for decoherence. It 
also determines the stochastic properties in the Langevin 
description of the stochastic limit. 

When the initial environment is a general Gaussian state, 
which may be generated by any combination of squeeze, 
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displacement, and rotation operators acting on the vac- 
uum state, the influence functional is, even for the ca se o f 
time-dependent coefficients, still given by the form (B2), 
with the noise and dissipation kernels expressed in terms 
of the Bogolubov coefficients which determine the time- 
dependence of the field operators when coupled to system 
variables. This fact allows the treatment of time-dependent 
background fields (e.g. electric fields) or curved spacetime 
problems that lead to an effective time dependence of the 
parameters in the Lagrangian. This explains the striking 
analogies between quantum optics and cosmological parti- 
cle productions p8fl . Hu and Matacz |34| derived a gener- 
alized influence functional and applied it to problems with 
arbitrary, but generally fixed, time dependence. Here we 
review the case where there are no cross-terms of the type 
<p a <p a , and the initial state of the field is the squeezed state 
given by 



Psqueezed (U] Ot) — S n (V (a) , (f) (a)) 
n 

X Pthermal 

(a) & n (r (a) ,</>(<*)), 



(B3) 



where pth is the initial thermal state for each environment 
mode: 

Pth(a) = ( 1 - cxp ( ) ) X! exp ' n ^~ Q 1 " " 



The squeeze operator is 
S(r (a), 0(a)) 

= iexp{r(a) (d'e"^ - a£e*««>) } , 



(B4) 
(B5) 



where a a , a J, are the annihilation and creation operators 
for each field mode (and polarization). 

The influence functional then has the form flB2|), where 
the dissipation and noise kernels are given by 



(B6) 



{ [a u (t) +&,(*)]* MO + MO] 
k (*)+M*)][Ou> (0+A* (t')]*} 



and 



v(t,t') 
1 



(B7) 



j , \ , / fko (ti 
2 J dwI ^ COth {-2hT) X 

{ cosh2r (lo) [o w (t) + & u (*)]* [00,(0 + ft, (*')] 
+ cosh2r (oj) [a w (t) + /3 U (t)] [a w (t') + /3 U {t')]* 
- sinh2r (lu) e' 2 ^ [at» (t) + /9 W (t)]* 
x MO + ft, (<')]* - sinh2r («) e 2i ^") 

x [oo, [o« (0+^(01 }■ 



The q;^ (t) , /3 W (£) are the Bogolubov coefficients determin- 
ing the time-dependence of the field annihilation and cre- 
ation operators. There are found by the solutions to the 
first order coupled equations 

a (t) = -iA (<)* 13 (t) - iB (t) a (t) , 

(5 (t) = iB (t) P (t) + iA (t) a (t) , (B8) 



where 



A(t) 
B(t) 



1 ( m a (t) u>l (t) m a {U) uj a (ti) 



2 \m a (ti) uj a (U) m a (t) 
1 ( m a (U) u> a (U) ( m Q (t) uj 2 a (t) 
2 



(B9) 



(*) 



m a (ti)u) Q (ti) 



APPENDIX C: CLOSED-TIME-PATH EFFECTIVE 
ACTION 

1. "In-in" generating functional 

In this appendix, we review the CTP method and use it 
to find the "in-in" one-particle-irreducible (1PI) effective 
action, from which we derive the real and causal equations 
of motion for the mean field. One may also use the CTP 
method to derive the "in-in" N particle irreducible (NP1) 
effective action, (N — > oo is called the master effective ac- 
tion), from which the Dyson-Schwinger equations of motion 
for the correlation hierarchy may be derived [^J?]] . This is a 
powerful tool for addressing problems where nonperturba- 
tive effects from higher order correlations need be treated 
in a self-consistent manner, and it may be applied to the 
problems we address here if dissipation and fluctuations are 
to significantly affect the mean-trajectory solutions. 

To formulate the CTP method, we define a "doubled" 
spacetime manifold M C tp = Mx{f,2}(46). Mis 
the Minkowski (or curved) spacetime manifold in between 
the initial time hypersurface E(tj) and the final time hy- 
persurface S(i/), where the two copies of M are joined. 
Spacetime points (x a ) carry a CTP index a, 6,... indicat- 
ing whether they live on M 1 or M 2 .At E/, where M 1 is 
joined to M 2 , the points (x 1 ,^) and (x 2 ,t 2 ) are identi- 
fied (see Figure 3). We will usually write (p a (x) instead 
of (p(x a ). M 1 (M 2 ) is often called the positive (negative) 
time-branch because positive frequency modes living in 
M 1 (M 2 ) evolve forward (backward) in time. The direc- 
tion of time of positive frequency mode propagation is 
determined by the sign in the exponential in ( pl[ ). One 
can keep track of this sign by introducing a CTP metric 
c ab _ §ai S bi _ §a2 S b2. c ab will a i ways used to contract 

CTP indices (e.g. (p a h a = (p x h\ — (p 2 hi). We define the 
CTP time-ordering operator T so that operators on M 1 
are time-ordered, operators on M 2 are anti-time-ordered, 
and all operators on M 2 are ordered left of operators on 
M 1 . Histories are to be thought of as beginning on M 1 at 
ti, moving forward in time to tf where M 1 joins M 2 at 
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£(£/), and then going backward in time to t\ on M 2 , thus 
the name 'closed-time-path'. 

We first consider a single free field. The "in-in" generat- 
ing functional is given by 

Z- m -in[h a ] = J dip f {<Pi\Vf) h2 (<Pf\<Pi) hl (Cf) 



= Tr 



Texp ' h ' dx ^H( x ) K(x) } Pi 



The subscript H indicates that the operators are in the 
Heisenberg picture in this expression. It follows that the 
CTP time-ordered expectation values are given by 



{T<p a H {x x ) - ■ - b H {x n )) Ja 

= Tr {f$j(xi) ■ ■■<p h H (xn)pi}j a 

6 n Z[h a ] 



(C2) 



Sh a (xi) ■ ■ -5hb(x n ) 



The generating functional has the path integral represen- 
tation 



z[K] 



CTP I n 



ip a h a ) pM) 



(C3) 



In the last line we used an abbreviated notation where 
S v [<p a ] = S^V 1 ] - S v [tp 2 ] , Pi (<p?) = and 
(p a h a = J dxtp a (x)h a (x). The CTP subscript on the inte- 
gral implies the CTP boundary conditions on integrations 
over the initial density matrix. We will assume that the 
initial state is approximated by the Gaussian 



Pi{9i) = -4exp 



^{x)K ab {x;x')^{x') 



where the quadratic kernel 



K ab (x 1 \x 2 ) = Kabiti,* 1 ;^,* 2 



(C4) 



(C5) 



vanishes everywhere except on £(ii). Any specific initial 
condition makes the evolution special, as it carries the par- 
ticular information at £j. It is in this restricted sense that 
Lorentz invariance of the theory is broken (unless p is the 
Lorentz invariant vacuum state for the full theory). Any 
non- vacuum state picks out a reference point, as is the case 
for the influence functional constructed with a specified ini- 
tial condition of the system and environment. There are 
many physical situations where the information of a chosen 
initial state or a fiducial time is encoded in the dynamics, 
and affects the attributes of a system. Most nonequilib- 
rium processes are of such a nature, while chaotic systems 
have particular sensitivity. For instance, if (C4) is a ther- 



mal state defined on the thermal fluctuations will 

only be isotropic in the rest frame determined by the time 
coordinate t. Another example which picks out a particu- 
lar reference frame is when the system and field environ- 
ment are initially uncorrelated. For interacting theories, 



exact factorizability holds only on a single spacelike sur- 
face since interactions inevitably lead to correlations both 
before and after such a special surface. For this reason, 
field-environment induced noise in a (particle) subsystem 
is not Lorentz invariant when one assumes that systems A 
and B factorize at the initial timep| 

We now consider the case where the field action S v [<p] 
has the quadratic form 



1 



S; [<p] = - / dx (d^ a d^p a + m 2 p aV a ) 



^ / dxdx'<f a (x)A ab (x,x')tp b (x') 



(C6) 



B.T., 



where 

A ab (x, x 1 ) = c a b {-d 2 + m 2 ) S (x, x) 

and B.T. is the b oun dary term. Further assumin g a Gaus 
sian initial state ( |C4|) , the functional integral in 
be solved exactly giving 

,-1/2 



Z[h a ] = (det (Gab))' 



) may 
(C7) 



x exp 



2h 



dxdx / h a (x)G ab (x,x')h b (x' 



where the Green's functions 



G ab (x,x') =Aj(x,x') 

S 2 S&a] 



(C8) 



S(p a (x)6(p b (x') 



are also given by 

G 11 G 12 
G 21 G 22 



-i (T + ip x ip X ') -i (tp x tp x >) 
-i(<Px'0x) -i{T~<p x 0x>) 

G F (x,x') G + (x,x') 
gJ[x,x') Gp(x,x') 



(C9) 



Gf is the Feynman and G F is the anti-Feynman (also 
known as the Dyson) Green's function and G^ are 
the positive (or negative) frequency Wightman functions. 
T + denotes time-ordering; T~ denotes anti-time order- 
ing. In the equations of motion the retarded (G ), ad- 
vanced (G A ) , and Hadamard (G H ) Green's functions play 
a prominent role. They are given by 



Gr(x,x') 
G a (x,x') 

Gh{x,x') 



-i6(x,x')([4> x ,^ x > 
-i6(x',x) (\4> x ,4> a 

-i ({^xAx'}) = G F 



Gf ~ 
= G F 

^Gp 



G+ (CIO) 
-G- 

— Gj- + G— . 



11 This observation is more generally true, of course. Only 
the dressed vacuum state of the interacting theory is Lorentz 
invariant, and in that case the influence functional cannot easily 
be found. So generically, the noise in the stochastic action will 
not be Lorentz invariant. Our belaboring on this general point 
arose from Prof. Ted Jacobson's insistence on clarity to his 
queries. 
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These various Green's functions are not independent. They 
satisfy the relation 



Gf ~\~ Gi 



G_. 



(Cll) 



Defining sum and difference variables h + = (h 1 + 
h 2 )/2, h~ ee (h 1 ~ h 2 ) and p+ = (cp 1 + p 2 )/2, ip~ = 
(p 1 — p 2 ), we can express the generating functional in terms 
of retarded and Hadamard Green's functions as 

= (detG ah )~ 1/2 exp{ - ± J dxdx' 

x [2h-(x)G R (x,x')h+{x') 

+ ih-(x)G H (x,x')h-(x')]}. (C12) 

The close relationship between the generating functional 
Z [h^] and the influence functional F [j^] is now apparent; 
in fact, 



^[^(detG"") 172 ^] 



(C13) 



when the source h is identified with the physical current j 
produced by the particles (r) . 

The generating functional for normalized expectation 
values, or connected correlation functions in the language 
of Feynman diagrams, is 



W[h a ] = -ih\nZ[h a ] 



(C14) 



-~ / dxdx'h a (x)G ab (x,x')h b (x') 



ih 



lndet G 



ab 



or in terms of the sum and difference variables h 
W[h a ] = - j dxdx'{h-(x)G R (x,x')h+(x') 



ih . 



+ l -h-{x)G H {x,x l )h-{x')\ - y lndetG afc 



== S 



ih 



IF 



•lndet G ab 



(C15) 



Thus, W [h a ] is just the Feynman- Vernon influence ac- 
tion, up to the log-det term. But for the quadratic action 
and initial state we are considering, G ab is independent 
of h, and so the log-det term is just a c onsta nt. Noting 
that h a (p a — h + <p- + h~ip + , we may use (C12) to find the 
mean-field and symmetrized Green's function as 



( V (x)) h = (p + (x)) h _ =0 = 



SW[h+,h- 



5h~{x) 



(C16) 



and 

(W(x),^(x')}) h = 



dxG R {x, x')h(x') 
S 2 W[h+,h- 



Sh (x)5h (x 



L = G H (x,x'). (C17) 



We will follow the standard terminology and call (p a [h a ] the 
classical fields, defined by 



<Pa[h a ] 



5W[h a ] 
Sh a 



(C18) 



This terminology is formal, since the physical (real- 
situation) meaning of classicality involves decoherence 
which requires coarse-graining. 

The one-particle-irreducible effective action is the Leg- 
endre transformation of the source dependent functional 
W^[/i a ] with respect to the classical field: 



T[$ a (x)] = W[h a ] 



dxh a (x)ip a (x). 



(C19) 



The real and causal equations of motion for ip a are found 
from 



ST[ip a ] 
5ip a (x) 



= -h a (x). 



(C20) 



For the quadratic case we are treating it is simple to invert 
the equations of motion giving 



h a {x) = / dx'G^{x,x')q> b {x') 
dx'A ab (x,x')tp b (x'). 



(C21) 



Substituting this result for h a (x) into (C19) gives 



Wa(x)} = i / dxdx'<p a (x')A ab (x,x')<p b (x') 



ih 



IndetG 



ab 



(C22) 



ih , 



Sim- ^IndetG 



ab 



Hence, we see that the effective action is just the classical 
action evaluated in terms of the solutions Cp a , plus a con- 
stant. This result is only the case for quadratic theories 
for which quantum corrections don't change the mean-field 
equations of motion. 

We can also proceed more formally, and generalize these 
results to non- quad ratic ( nonli near) field actions. First, 
we substitute ( C2C| ) into ( |C19| ). Using the path integral 
representation (CE) for Z[h a ], we find the functional inte- 
grodiffcrential equation for the effective action 



ihlni D(p a exp-(S[ip a ] 



'CTP 



r [tpa] = -ihln\ j 
Jc 

- J dx—(<p a - (p a )) Pl {ip a )\ 



(C23) 



Defining the fluctuation field (p a {x) = <p a {x) — <p a {x), and 
using Dip a = Dip ai we may rewrite this as the functional 
integral of the fluctuations around the classical fields (p a : 



I CTP 



D<p a exp-(S[<p a + p a ] 



T [(f a ] = -ih\n\ 

Jc 

- J dx-^-p a )Pi{(Pa + (f>a)\ 



(C24) 
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Equation (C24) has the formal solution 



2. CTP coarse-grained effective action 



T[(p a (x)} = S[<p a ] - ylndet (G ab (x,x')) + itiT 1 [<p a } 



where 



G ab (x, x')(p — 



S 2 S* [if] 



6lp a (x) 6lfi b (x ; ) 



(C25) 



(C26) 



<Pa = Va 



The G ab (x, x')ip are the one- loop propagators for the fluc- 
tuation fields (p a , but now they are functionals of the clas- 
sical fields <p a . The functional Ti is given by the sum of 
one-particle-irreducible (1PI) graphs with external legs de- 
termined by pi (</?), and with propagators given by the one- 
loop fluctuation field propagators. These quantum correc- 
tions begin at two-loop. Vertices are given by the shifted 
interaction term for the fluctuation fields dp : 



Sint[Pa\ = S[ip a + Va] ~ S[lf a ] 



SS[<p a ] . 
dx — ip a (C27) 

0<Pa 



j dxdx' 5 S ^ a \ (p a {x)(p b {x'). 

The one-loop contribution to the effective action is given by 
the log-determinant term, but now, unlike for the quadratic 
action case, this term is not independent of ip a . For one- 
loop calculations one neglects F%, however, under some cir- 
cumstance important nonperturbative (with respect to the 
coupling constant) effects only begin at two- loop ps[| . For 
instance, the time-dependent Hartree-Fock approximation 
involves two-loop corrections. 

Both the influence functional and CTP approaches may 
be used to find the full quantum theory of the particles, and 
not just the semiclassical limit. Within the CTP frame- 
work, this may be done in the context of the N-particle- 
irreducible effective action, or its full extension (N — > oo), 
the master effective action The full quantum theory is 
equivalent in content to the Dyson-Schwinger hierarchy of 
correlation functions Q . By using the master effective ac- 
tion, this hierarchy may be truncated consistently at finite 
order by slaving higher order correlation functions to some 
finite number of lower order correlation functions. Within 
the influence functional approach, one may derive the finite 
time propagator for the system's (in these studies, parti- 
cle's) reduced density matrix. This propagator then pro- 
vides a description of the effective dynamics of the particles 
with backreaction from quantum fields. This may then be 
used directly, or as an intermediate step, in deriving the 
master equation. Solving for the prop agator requires eval- 
uating the functional integrals in ( |2.8[ ). For linear theories, 
this may be done exactly. For nonlinear theories, one may 
either expand the functional integrals around the station- 
ary phase solution giving the loop expansion (in order of 
H) for the propagator, or one may exp and the influence 
functional in the coupling strength ( |2.9| ), giving coupling 
constant perturbation theory for the reduced density ma- 
trix propagator. 



The coarse-grained effective action (CGEA) was first in- 
troduced by Hu and Zhang |T^] to treat the backreaction 
of one subsystem on another, where their separation is jus- 
tified by the existence of a physical parameter (e.g., heavy- 
light mass, fast-slow time, high-low modes) which one uses 
to carry out a saddle-point expansion for the action of the 
open system. It includes the ordinary (n-loop) effective 
action as a special case where the expansion parameter 
is h. Here we apply CGEA to interacting particles and 
fieldspl In this sub-appendix, we extend the Closed-Time- 
Path method to the particle plus field model. 

To produce a generating functional for both particle and 
field correlation functions we couple both particle and field 
variables to independent sources, adding the term 



E 



dTJ nf j, (r) z% (r) + / dxh (x) ip (x) 



(C28) 



to the action. For N particles (indexed by n) in D space- 
time dimensions (indexed by /i) we must have NxD sources 
J nAt independently coupled to the spacetime coordinates of 
each particle. However, rather than considering the full 
generating functional for all particle and field correlation 
functions, we shall study the open system dynamics of 
particles with the influence functional obtained by coarse- 
graining (integrating out) the field degrees of freedom. We 
do this by not including a field-source h (x) . 

The exact "in-in" connected generating functional (for 
particle correlation functions) is then given by 



W[J a ] = -ih\n[ [ D<p a Dz a exp{?-S a {<p a , z a ] 



CTP 



(C29) 



where 



s[<p a ,z a ] = s;[<p] + s a [z] + s? nt {<p,z} 

Sfnt [<P, A = J dxf (x) ip a (x) 

S a [z] =S[z]-S [z'} , (C30) 



and S" [tp] is the quadratic field action given in 



We 



are considering cases like QED, where the nonlinearity arise 
only through interactions between the particles z and field 
ip. As per our earlier discussion (see Appendix A), the path 
integrals over the worldlincs arc restricted to U < z°. The 
CTP temporally-ordered correlation functions (Tz-iZ2---Z n ) 



This Appendix is formal, in that we assume a worldline 
restricted-path integral structure without developing the full de- 
tails, which is the topic of our second series. 
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are found from the functional derivatives of W[J a ] with 
respect to J a (r) . We define the "classical" particle solutions 
to be 



8W[J a ] 
8J a {r) 



(C31) 



The solutions z a [J a ; y) are functionals of the sources J a . As 
we mentioned earlier, calling z a the classical solutions is a 
formal terminology since classicality involves other issues 
like decoherence. The CTP coarse-grained effective action 
is defined by the Legendre transform 



r[*«(T)] = W[J a ] - / dTj a (r)z a (T). 



(C32) 



r[z a (r)] is the coarse-grained one-particle-irreducible gen- 
erating functional for particle correlation functions. We 
find the real and causal equations of motion for z a from 



ST[z a ] 
Sz a (r) 



J a (r). 



(C33) 



When J 1 = J 2 , the solutions on the two time- branches 
M a are equal. The physical solutions are found by setting 
J a = 0. 

We may now use the results from the previous Appendix 
to integrate out the field tp a keeping z a fixed, and assuming 
that the initial field and particle states are factorizable. We 
simply replace the source h a (x) with the physical current 
j a (x) . The connected generating functional is then given 

by 

W[J a ]=~ih\n[ / Dz a exp^-{s a [z} + S IF [j a ] 
Jctp n, 

: Tr In G ab + £ / drJ^z^} p(zf)} . (C34) 



f , 



The trace-log term just comes from the normalization fac- 
tor (detG ab ) wh ich is inde pendent of z. 



Substituting (C32) into ( |C32| ) gives an integrodifferential 
equation for the effective action 



T[z a ] = -ihln [ f Dz a e^{s a [z] + S IF [j a ] 
Jctp " 



-TrlnG Qb 
2 



(C35) 



jdr§- a {z a -z a )} P {zt) 



The difference z a — z a is the deviation of the particles path 
integration histories from the classical trajectory z a . Shift- 
ing variables by defining the fluctuation coordinate 



z a (r) = z a (r)-z a (r), 



(C36) 



and using the invariance of the functional m easur e under 
translation: Dz a — Dz a , we may express (C35) as the 



functional integral of fluctuations z a around z a . We obtain 
T[z a ] = ~th\n [ f Dz a exv-{s a CGEA [z a + z a ] 

- iTrlnG ah - J dr ^-J*} p{z«)] , (C37) 



where S^ GEA [z a + z a ] = 5° [z a ] + S IF [z a + z a ] is the coarse- 
grained effective action written in terms of the new vari- 
ables z and z. 

Now the restriction on the particle-worldline path inte- 
grals means that even the one-loop approximation goes be- 
yond the exact Gaussian integrations. This is a finite-size 
effect, which we discussed earlier. So that we can avoid 
considering finite-size corrections here, we will assume that 
the initial particle state is defined in the infinite past, and 
hence the path integrations over the time coordinate z° (t) 
have the full range from minus to plus infinity. Then the 
effective action has the formal solution 

T [z a ] = Scgea [z a ] - y lndet G ab - y lndet (B^ 1 ) 



where 



iB ab (T, t') 



>CGEA 



5z a {T)8z b (r') 



(C38) 



(C39) 



is the one-loop propagator for the particle fluctuations z a 
around the classical solutions. Scgea [z a ] is the same ef- 
fective action found from the influence functional method. 
Ti [z a ] contains two-loop and higher quantum particle cor- 
rections that arise in nonlinear theories. It is given by 
the sum of all one-particle-irreducible graphs with exter- 
nal legs determined by pi(z a ), internal propagators given 
by B a i,(T, t')z, and vertices given by the shifted interaction 
term for the fluctuation coordinate z. 



S ln t[z a ] = S C GEA[z a + Z a ] - S CG EA[z a 
( SScGEA r _ Q A xa 



d:T 



\ Sz 



^ [z a ] z a (r) 



-- [dr [dr'( S2 , S °? E h ? ; [z° 
2 J J \8z a (T)5z b (T>) [ 

X Z a ( T )z b {T'). 



(C40) 



From (C39) and (C40), it follows that the sum of graphs 
for Ti have vertices and propagators that depend on z a . In 
terms of Feynman diagrams, the gr aphs that are included in 
Ti have only particle propagator (C39) lines. The implicit 
dependence in the graphs on the integrated field variables 
comes about through the dependence of the vertices on the 
solutions z a , which in turn depend implicitly o n the average 
field properties through z's defining equation (C33). 
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APPENDIX D: CORRESPONDENCE BETWEEN 
LANGEVIN AND HEISENBERG EQUATIONS OF 
MOTION 

The construction of the stochastic effective action made 
use of a Gaussian identity to rewrite the real part of the 
influence functional as the stochastic average of the char- 
acteristic function of a noise. But we might still ask if all 
aspects of the field's quantum fluctuations are encoded in 
the statistics of the noise, or if different states of the quan- 
tum field can give the same noise in the stochastic limit. 
This possibility has been pointed out by Raval in |m| . For 
a Brownian particle linearly coupled to a scalar field, the 
Heisenberg equations of motion are 



dt 2 



q(t) + LO 2 q(t) 2 = e0(x(t)) 



and 



0%v>(x(t),t) = -$(t)a(x-x(t)), 



(Dl) 



(D2) 



where the particle, with internal degree of freedom q, mo ves 
on a fixed trajectory x(i). The causal solution to (D2) in 
terms of the retarded Green's function plus a homogeneous 
free field solution gives 



d 2 „, , 
dT 2 ^ 



dt'G R (t,t')q(t')+LU 2 q(t) = o (x(t),t). 



(D3) 



When G R (t,t') is local (as is the case for massless fields), 
we may integrate by parts and write 



d 2 d 

^2 $(*) + ie J t ^ + W 0$(*) = £o(x(t), t) 

with the solution 

q(w) = (ijJ 2 + ieuj - ujq) 1 ipo(uj, x(t)) + qo(ui) 
= 9(u)0o(u,x(t)) + q {uj), 

where 

0o(w,x(t)) = / *e~Vo(x(t),t). 



(D4) 



(D5) 



(D6) 



We may now use (D5) to find both the particle commuta- 
tors (neglecting the homogeneous solutions qo(u))), 

([q(u;'),q(cj)})=g(cj')g(u;) (D7) 
x ([^^t)),^,^))]) 

and the particle anticommutators, 

(D8) 

x ({£(£,/, x(t)),0(a;,x(t))}). 

The corresponding Langevin equation rep laces the homo- 
geneous field operator solution <po in (D3) with the noise 



£(s), and the quantum operators q(t) are replaced with the 
stochastic variables q$(t). The solution otherwise follows 
as above, with 



9{M = g(w)£(a;,x(t)). 



(D9) 



Again, we drop the homogeneous solution to q^(u>). The 
symmetrized commutator of q^(u>) is then given by 

({ ?e (,'),^H})^( w ')sH 

x ({£(<■/, x(t)),£(a;,x(t))}) 

= g{u')g{w) (D10) 

x ({0(c./,x(t)),0(w,x(i))}), 

where in going from the second to the last step we have used 
the equality of the symmetrized noise and field operator 
commutators. But 

= SM (K(^',x(t)),^,x(t))]) 
= (Dll) 

^([ g V),?M])- 



Equation (Dll) follows as a consequence of 



(D12) 



Therefore, it appears to be formally possible for the Heisen- 
berg equations of motion to differ from the Langevin equa- 
tions by terms proportional to the field commutator. The 
origin of this possible discrepancy between the stochas- 
tic and Heisenberg correlation functions remains unclear. 
Raval has speculated that it may be related in some way 
to the extremization principle used to find the stochastic 
effective action Q|. Because (<£o) = (0 = 0, the mean 
equations of motion will exactly agree, so any difference 
between the stochastic and quantum theory will appear in 
the higher order fluctuations. 
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Particle/Field Open Systems 




FIG. 1. A schematic depiction of the relationship between 
the classical, semiclassical, stochastic, and quantum regimes. 




FIG. 2. Example of two oppositely charged particles in a 
background electric field, whose semiclassical motions z e ± are in 
causally disjoint regions of spacetime. The radiation force R12 
vanishes, but the correlation-noise Nn does not (for explanation 
see Sec. 5). 




FIG. 3. The closed-time-path manifold Mctp is composed 
of two copies of the ordinary spacetime manifold, Mi and M2, 
joined at a final-time hypersurface E (tf) . 
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